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Pinal  Scientific  Report 


A  conference  summarizing  ten  years  of  work  on  the  Singularity 
Expansion  Method  (SEM)  and  the  Eigenmode  Expansion  Method  (EEM)  was 
held  at  the  Carnahan  House  of  the  University  of  Kentucky  in  November 
1980.  The  proceedings  of  this  conference,  originally  due  for 

publication  in  Fall  1981,  did  not  appear  until  the  Spring  of  1982 
in  Volume  1,  Number  4  in  "Electromagnetics"  (October-December 
1981)  . 


These  proceedings  contain  two  papers,  one  by  C.  L.  Dolph  and 
one  by  A.  G.  Ramm  which  reviews  work  on  this  subject.  These  papers 
are  reproduced  in  Appendix  I  and  II  respectively  of  this  report. 


Subsequently  C.  L.  Dolph  submitted  an  expanded  version  of 
A.  G.  Ramm's  paper  to  the  Journal  of  Mathematical  Analysis  and  Ap¬ 
plications  [Vol.  86,  No.  2,  April  1982],  This  paper  was  entitled 
"Mathematical  Foundations  of  the  Singularity  and  Eigenmode  Expansion 
Methods  (SEM  and  EEM).  In  this  author's  opinion  it  represents  the 
best  available  mathematical  treatment  of  the  subjects  it  covers.  It 
is  reproduced  here  as  Appendix  III. 


Section  5  of  this  paper  entitled  Problems  is  still  viable. 
There  have  been  claims  that,  under  some  circumstances,  the  complex 
roots  of  the  Green’s  function  for  the  exterior  Dirichlet  or  the 
Neumann  Laplacian  are  simple.  To  this  author's  knowledge  this  is 
still  an  open  problem. 
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Although  Ramin  and  Dolph  have  repeatedly  emphasized  that  the 
Picard  method  is  only  valid  for  a  restricted  class  of  surfaces,  all 
the  work  that  this  author  has  seen  on  equivalent  circuits  is  based 
on  the  validity  of  the  Picard  method  and  consequently  must  be  con¬ 
sidered  of  limited  value. 

Two  other  topics  have  been  treated  since  the  above  conference. 

The  first  of  these  appears  in  a  paper  "Converence  of  the  T-matrix 
Approach  to  Scattering  Theory"  ,  Journal  of  Mathematical  Physics 
Vol.  23  (6) ,  June  1982,  and  its  subsequent  generalization  in  a 
preprint  by  the  same  title,  coauthored  by  G.  Kristanson,  A.  G.  Ramm, 
and  S.  Strflm.  This  represents  work  done  at  the  Institute  of  Theo¬ 
retical  Physics  in  Goteborg,  Sweden,  in  the  summer  of  1982.  In 
contrast  to  the  first  paper  on  this  subject,  the  second  paper  was 
not  done  under  the  auspices  of  this  grant. 

These  papers  are  reproduced  as  Appendices  IV  and  V,  respectively. 

One  large  area  still  under  investigation  involves  the  use  of 
variational  principles  in  these  problems.  It  will  perhaps  be  re¬ 
called  that  M.  S.  Agranovitch  carefully  avoids  these  even  though 
he  laid  much  of  the  foundations  for  the  EEM  and  SEM  in  his  appendix 
"Spectral  Properties  of  the  Diffraction  Problems"  which  is  contained 
in  the  book  Generalized  Methods  of  Normal  Modes  in  Diffraction 
Theory  by  N.  Voytovich,  B.  Katsenelenbaum  and  A.  Sivov,  Moscow 
(1977) .  More  specifically  he  states  "the  formal  Ritz  method  of 
finding  stationary  values  of  the  functions  (see  Chapter  III)  was  not 
analyzed"  . 
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V 


This  problem  has  been  around  for  a  long  time.  A.  G.  Ramm,  in 
"Variational  Principles  for  Resonances,  II"  ,  Journal  of  Mathema¬ 
tical  Physics  Vol.  23  (6),  1982,  developed  an  interesting  approach 
similar  to  that  used  in  his  paper  on  the  T-matrix.  This  paper  is 
reproduced  as  Appendix  VI  . 

This  author  has  been  exploring  min-max  theory  in  the  hope  of 
obtaining  error  bounds  for  scattering  problems.  Earlier  papers 

"A  Saddle  Point  Characterization  of  the  Schwinger 
Stationary  Points  in  Exterior  Scattering  Problems", 

J.  Soc.  Indust.  Appl .  Math.  Vol  5,  No.  3,  September  1957. 


"The  Schwinger  Variational  Principles  for  One-Dimensional 
Quantum  Scattering"  (with  R.  K.  Ritt) ,  Math.  Zeitschrift 
Band  65  (1956) ,  309-326. 

If  the  orientation  of  the  saddle  could  be  determined  a  priori,  it 
appears  possible  to  use  relative  cycle  theory  to  obtain  estimates. 
A  search  of  the  literature  has  not  been  helpful  in  that,  while 
min-max  theory  is  highly  developed  in  the  theory  of  games  and  in 
control  theory,  fuzzy  set  theory,  and  critical  point  theory,  no 
error  estimates  appear  to  be  known. 

An  earlier  paper 

"Symmetric  Linear  Transformations  and  Complex  Quadratic 
Forms"  (with  J.  E.  McLaughlin  and  I.  Marx),  Comm.  Pure 
Appl.  Math.  Vol.  VII  (1954),  621-632 
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still  has  not  been  generalized  to  the  infinite  case. 


Two  other  earlier  papers 


"A  Critique  of  Singularity  Expansion  and  Eigenmode  Expansion 
Methods"  (with  V.  Komkov  and  R.  A.  Scott) ,  Proc.  of  Confe¬ 
rence  on  Acoustics ,  Electromagnetic  and  Elastic  Wave- Scattering 
Focus  on  T-Matrix  Approach,  Pergamon  Press  (1980),  453-462. 


"On  the  Relationship  Between  the  Singularity  Expansion  Method 
the  Mathematical  Theory  of  Scattering"  (with  S.  K.  Cho) ,  IEEE 
Trans,  on  Antennas  and  Propagation,  V.  AP-28,  No.  6  (1980), 
888-897. 


are  still  pertinent  to  any  evaluation  of  the  SEM  and  EEM  methods. 
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APPENDIX  I 


ON  SOME  MATHEMATICAL  ASPECTS  OF  SEM,  EEM 
AND  SCATTERING 


C.  L.  Dolph,  The  University  of  Michigan.  Department  of  Mathematics. 
Ann  Arbor,  Ml  48109 


ABSTRACT 


The  relationship  between  the  integral  equations  usually  used  in  SEM  and 
the  scattering  matrix  is  examined.  Alternate  integral  equations  which  exhi¬ 
bit  only  the  poles  of  the  S  matrix  are  given.  Examples  are  used  for 
illustration  for  a  solvable  case. 

The  analytic  Fredholm  theorem  in  Banach  spaces  is  discussed  and  its 
advantages  for  numerical  calculations  emphasized. 

The  relationship  between  EEM,  SEM  and  the  theory  of  nonselfadjoint 
operators  is  briefly  discussed. 

INTRODUCTION 


The  ideas  lying  behind  the  Eigenmode  Expansion  Method  (EEM)  appear  to 
have  been  introduced  for  the  first  time  by  Kacenelenbaum  in  1969  [6.4] . 

The  Singularity  Expansion  Method  (SEM)  was  first  introduced  by  Baum  in 
1971  [3.1]  and  shortly  thereafter  independently  he  introduced  EEM  [3.4] .  The 
best  review  paper  of  these  USSR  contributions  is  that  due  to  Voitovic, 
Kacenelenbaum  and  Sivov  [6.11]  and  that  of  the  USA's  contributions  (in  this 
author's  opinion)  is  that  of  Baum  [2.2).  The  most  complete  review  of  the 
Russian  work  through  1976  is  the  Russian  book  [6.11]  by  the  above  three 
Russian  authors.  This  book  also  contains  a  mathematics  appendix  by 
M.S.  Agranovic. 

A  glance  at  the  official  bibliography  makes  it  clear  that  extensive 
work  has  been  undertaken  and  completed  since  these  beginnings,  and  more  will 
be  discussed  in  these  preceeding. 

In  view  of  the  extensive  publications  the  author  thought  it  might  be 
most  useful  to  provide  a  brief  guide  to  some  of  the  recent  mathematical 
developments  without  excessive  detail  and  without  proofs. 


This  work  was  partially  supported  by  the  Air  Force  Grant  800204 

Elect  romagnttici  1:375-383, 1981 

0272-6343/81  /040375-09$2.25 

Copyright  e  1981  by  Hemisphere  Publishing  Corporation 
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SEN 

For  the  scaler  wave  equation  an  exterior  Dirichlet  problem  would  be 
formally  given  by  the  first  four  of  the  following  equations.  The  fifth 
equation  is  the  solution  given  in  terms  of  generalized  eigenfunctions  which 
are  distorted  plane  waves  playing  the  role  of  the  plane  waves  used  in  the 
Fourier  integrals  which  occur  when  no  obstacle  is  present.  The  functions 
<x(k )  and  SC-.!  are  related  to  the  initial  conditions.  This  last  formula 
has  been  rigorously  established  by  Shenk  [6.92}  in  a  manner  similar  to  that 
used  by  Xkcbe  [s.3sj  for  the  quantum  mechanical  case.  Explicitly  the  genera¬ 
lized  eigenfunctions  are  defined  by  (6),  (7)  and  (8).  As  will  be  discussed 
below  several  different  methods  are  available  for  the  construction  of  V 
,2 

Cl)  AV  »  ~ 

at 

(2)  u(x,0)  -  ix 

(3)  |H  (x,  0)  =  f2 

(4)  U  *  0  on  r 


(5)  C(x,t)  =  - ~- 

(2c)J/^ 

(6)  (A  +  k3)i*  -  0 


(7}  -  0  on  r 

(8)  «+  -  elk-'  +  Vf  <x,k) 

In  contrast  to  the  operator  theory  approach  employing  the  continuous 
spectrum  SE>1  employs  the  Laplace  transform  which,  after  a  suitable  rotation 
in  the  s-plane,  can  be  defined  by  (9),  (10),  (11),  and  (12).  Condition  (10) 
is  one  form  of  the  radiation  condition  which  is  need  to  guarantee  uniqueness 
of  all  It  ,  I»  k  >  0  . 

(9)  V(x,k)  £  [  U (x, t)  eikt  dt 

•0 

(10)  AV  +  k2  V  =  -£ 

(11)  V  *  0  on  r 

(12)  — -  -  ikV  <=  0(|x|_1)  . 

3lx! 

The  function  V(x,k)  is  sought  in  terms  of  the  Green's  function  as 
given  in  equation  (13).  The  Green’s  function  in  this  equation  is  not  the 
well-known  Green's  function  of  free  space  but  is  determined  by  (14),  (15), 
(16),  and  its  domains  of  analytic  and  mernmc rphicity  are  given  by  the  next 
two  statments  (17)  and  (18).  f  o  Dolph,  M'-'  od  anil  Thoe  (6.25). 


^(x,k)(a(k)  eikt  +  8  (k)  e* 


ikt. 


d3k 
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(13)  V(x,k)  =  |  G(x,y,k)f  (y)dy 


where  G  satisfies 


(14) 

(15) 

(16) 

(17) 

(18) 


(&  +  k  )G  =  -6  (x-y)  in  0 


G  *  0  on  r 
3G 


ikG  =  Odxp1) 

3|x| 

G(x,y,k)  analytic  In  k  >  0 
G(x,y,k)  meromorphic  Im  k  <  0 


Once  V  (y ,  k)  has  been  found  the  solution  of  the  original  problem  can 
be  given  in  terns  of  the  inverse  Laplace  transform. 


If  a  >  0.5,  Im  k  >  -b,  b  >  0  and 

(19)  |  Vj  <  - - — -  ,  |  Re  kj  -*•  »  . 

1+  |k|a 

Then  for  0  <  Y  <  b 


i  riY*“  -ikt 

U(x,t)  =  V(x,k)e  1KC  dk 

J-iy-o 

Pushing  the  contour  down  yields  [6.69] 

n  -ik  .t  —  1 1  k.  1 1 

(20)  U(x,t)  «  l  e  3  V(x,k.)  +  0(e  "  n  ) 

j=l  3 

The  function  V  and  k's  which  occur  in  this  asymptotic  formula 
are  the  complex  eigenfunctions  and  eigenvalues: 


(21) 

(22) 

V. 

J 


6V.  +  k?V .  =  0  in  fi 
J  3  3 

=  0  on  r 

grows  exponentially  in  x  . 

Several  comments  are  now  in  order: 


(i)  The  estimate  (19)  is  valid  for  the  Dirichlet  problem  if  the 
body  is  (a)  star-shaped  and  (b)  non- trapping  in  the  sense  of  Lax  and 
Phillips  (6.52). 

(ii)  The  method  is  not  very  useful  since  it  involves  the  construction 
of  the  Green's  function  and  then  the 'determination  of  its  poles. 

(iii)  It  is  an  open  problem  to  find  conditions  when  the  asymptotic  series 
(20)  will  actually  converge. 


Instead  in  SEM  it  is  usual  to  employ  the  methods  of  potential  theory. 
For  the  exterior  time-independent  Dirichlet  problem  corresponding  to  the 
time  dependent  problem  we  have  been  considering  up  until  now,  this  involves 
consideration  of  the  following  set  of  equations  whichemploy  the  known  Free 
space  Green's  function. 
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(&  +  k2)V  =  0 

V  =  -V.  on  r 
me 

ikjx  -  v| 

*0<*’*>  =  h  T£~£P 

f  3 

(23)  V<x,k)  =  jTT  <>o(ii'il)d(il)d3v 

Jr  y  Y 

(24)  I  +  8(k)  d(v)  +  X  -  ^(x.yOd^Jd  y  =  0 

X  =  X (k)  =  -1 

Using  the  Fredholm  alternative  the  poles  are  sought  as  non-triviai 
solutions  of  the  homogeneous  integral  equation  (24) .  Those  which  may  occur 
for  real  k  correspond  to  eigenvalues  of  the  associated  interior  Neumann 
problem.  As  such,  as  we  shall  see,  they  occur  because  of  the  double-layer 
assumption  and  can  be  eliminated  by  other  assumptions.  As  shown  by  Dolph 
and  Wilcox,  see  Dolph  (6. 96)  they  do  not  contribute  to  the  scattered  field 
nor  do  they  appear  in  it  for  any  separable  case. 

The  homogeneous  integral  equation  which  occurs  here  can  be  treated 
mathematically  several  different  ways.  Marin  [3.9]  employed  Carleman’s 
Hilbert  space  theory  but  the  analytic  Fredholm  theorem  attributed  to 
Steinberg  [6.S4]  is  perhaps  the  most  convenient  since  it  is  applicable  in 
mors  general  Banach  spaces.  Since  matrix  approximations  are  used  in  the 
numerical  calculation  of  the  poles  the  choice  of  the  Banach  space  of  con¬ 
tinuous  functions  is  perhaps  the  most  convenient.  See  Dolph  and  Cho  [6.22] 
for  a  fuller  discussion.  For  a  Hilbert  space  the  proof  of  the  analytic 
Fredholm  theorem  can  be  found  in  Reed  and  Simon  [6.84). 


Analytic  Fredholm  Theorem.  -  Steinberg  [6.94] 


Let  0(B!  =  set  of  bounded  operators  on  the  Banach  space  and  let.  K  be 
an  open  connected  subset  of  the  complex  plane.  T (K)  is  analytic  in  K  if  for 
each  k.Q  €  K 


v  T  (k-k  >n 
0  n  0 


T  e  0(B) 
n 


Theorem.  If  T(k)  is  an  analytic  family  of  compact  operators  for  k  6  K  , 
then  either  1  -  TU)  is  nowhere  invertible  in  R  ,  or  else  [I-T(k))-*  is 
cveromorphic  in  K  . 


If  B  is  a  separable  Hilbert  space,  the  residues  are  finite  rank 
operators. 


One  way  of  eliminating  the  poles  which  are  not  intrinsic  to  the  exterior 
scattering  problem  is  to  replace  the  Ansatz  (23)  of  the  double  layer  by  the 
corip; -x  combination  of  a  double  and  single  layer  an  used  by  Brakhage  and 
Werner  for  the  Dirichlet  problem  (6.11)  and  by  Kussmaul  [6.50)  for  the 
Neumann  problem.  In  the  latter  case  additional  difficulties  need  to  be 
overcome  because  of  the  high  order  of  the  singularity. 

For  the  Dirichlet  Problem  the  Ansatz 
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“  f  ( —  iT  V*'*,u<z>dv° 

J  i  y  7 

leads  to  the  homogeneous  integral  equation: 

(25)  U(x)  +  |  (  - it)  *0<*ii)W  <X)dy°  =  0 


T(X)  =1 

for 

Re 

k  > 

0 

*=  0 

for 

Re 

k  < 

0  . 

This  equation  has  only  trivial  solution  for  I«  k  >  0  and  hence  the  only  non¬ 
trivial  solutions  can  occur  for  Im  k  <  0  and  as  Ramm  (6.72)  has  shown  these 
occur  at  the  poles  of  a  Green's  function  and  are  in  fact  the  intrinsic  poles 
of  SEM. 


The  non-trivial  solutions  of  this  last  equation  for  Im  k  <  0  also  agree 
with  the  poles  of  the  S  matrix.  The  S  matrix  is  generally  thought  to  con¬ 
tain  all  intrinsic  properties  and  in  fact  Lax  and  Phillips  have  given  two 
proofs  of  the  fact  that  the  S  matrix  uniquely  determines  the  obstacle  for 
the  Dirichlet  problem  —  see  Theorem  (5.6)  cr.  [6.52],  Chapter  V. 

For  the  problem  here  it  can  be  shown  that  the  V_(x,h)  of  (8)  and  the 
S  matrix  are  related  by  the  formulas:  [In  the  last  equation  the  integral 
operator  is  compact). 

X  =  r£,  5.  =  kw 

-ikr 

V_  (r9  , kw)  =  s—^—  (s_(9,k,w)  +  0(D) 

S  (k)ra(  )  =•  n(6)  +  ^  I  m(w)s  (6,k,w)*dSw 

2Mw|=i  ‘  .  . 

The  complex  eigenvalues  are  poles  of  S(k)  . 


Derivations  of  these  formulas  can  be  found  in  Lax  and  Phillips  [6.52) , 
in  Schmidt  [6.87]  for  the  quantum  case  of  the  Schrodinger  equation  and  for  a 
very  general  case  in  Shenk  and  Thoe  (6.9).  A  physical  derivation  of  the  last 
formula  is  due  to  Saxon  [6.86]  is  also  contained  in  Dolph  and  Cho  [6.22J. 

This  last  paper  also  contains  an  appendix  in  which  a  heuristic  derivation  of 
the  mathematical  theory  of  scattering  initiated  by  Jauch  is  given. 


For  the  cylinder  (24)  becomes 
ika 


"<a»V  +  4 


\Q  F(!l7  Ho1><ka  sin'#)  <a*Wd3o 


and  has  solutions  given  by 


i(a,«  )  -  l 


2(-l)n-1Jn(k  a)e 


inf. 


—  J1(ka)H,1) <ka) 
n  n 


0 


The  complex  roots  of  the  Hankol  function  are  intrinsic,  those  of  the  deriva¬ 
tive  of  the  Bessel  function  well-known  to  be  those  of  the  associated  interior 
Neumann  problem. 
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The  Brakhage- Werner  equation  corresponding  to  (25)  is 


V 


ika  f2T_l_ 

4  k  b<ka)  "  iXJ 


sin  6, 


ka 


°‘a'W 


0 


with  a  solution  exhibiting  only  intrinsic  poles;  namely 
,  ■»  (-l)n_1J  (ka)e 

O(a'-o)  '^T  l  ~r~) - I -  * 

"  -»  (ka)  (J  {ka)  -  iJ  (ka) ) 

n  n  n 

For  this  problem  the  cor.piex  eigenfunctions  are 
V(r,6 )  =  e1"-  Ha)  (k.r) 


where 


(A  +  kQ)V  *0,  V  -  0,  r  =  a 


and  the  scattering  tutrix  is  given,  as  shown  by  Shenk  and  Thoe  (6.911  to  be 
S(k)  [l  ay^)  =  -  l  -77 


_  u(2)  .  im9 

®  H  (ka)a  e 


lr.  most  c.-.st-s  it  is  necessary  to  resort  to  matrix  approximation  or  to 
have  methods  for  the  calculations  of  the  poles.  In  the  case  of  the  former, 
Rar.n  (6.72]  has  Established  the  following: 


Foies  coir.oi.de  with  fc 
is  not  invertible.  Let 

in  H  ---  L2  ( T )  .  Then  if 

.  n 


for  which  I  +  B(k)  of  (24) 


{f.> 

3 


be  an  orthonormal  basis 


y  c  .  f  . 


1  3  3 

<  (X*B(k))fi, 


f  .  > 
3 


It  follows  that 

n 

Y  b .  .  (k)  c  .  =  0 

jtl  13  3 

(n)  ■  •’* 

Let  k  ,  n  «*  1 , 2, 3. . . ,  be  the  roots  of 


det  b^ (k)  «  0 

Then  the  limits  Lite  k^  *  k  exist  and  arc-  the  poles  of 

the  ps.es  of  the  Greens  function  G  .  Every  pole  of  G  can  • 
be  obtain  in  this  way. 


BEM 

For  this  s-C"  time- independent  Dirichict  problem  the  Eigenmode  Expansion 
method  would  involve  the  following: 
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♦  k2  «  0 


A(k)$ 


g  on  T 

ik|  x-y| 


4n|x-yj 


$ (y)dy 


Ansatz 


A(k)$  «  g  on  T 

A(k)$  <=  X  <k)$  on  T 

n  n  n 


Picard  method  gives 

i.lllpL 

X  *n 

when  is  this  valid? 


The  Picard  process  is  certainly  valid  for  the  cylinder  and  the  sphere.  In 
fact  as  first  noted  by  Kacenelenbaum,  Sivov  and  Voitovic  [6.11]  for  the  cylin¬ 
der  they  are  explicitly  given  in  the  case  of  even  0  by 

$  (8)  ”  cos  n8 
n 

X  Ik)  =  ~  H  <ka)J  Cka) 
n  z  n  n 

$  (r,  )  =  H^(ka)J  (kr)cos  n9  ,  r  <  a 
n  ft  n  — 

(kr)  J  (ka)  cos  n9  ,  r  >  a 


While  Dolph  [6.20]  appears  to  have  been  the  first  to  suggest  the  use  t 
non-self-adjoint  operators  in  scattering  problem  Agranovic  in  [6.1],  [6.2], 

[6.3]  and  Raima  [6.72],  [6.73]  appear  to  be  the  first  to  systematically  apply 
this  idea.  Raram  in  particular  considered  the  Hilbert  space  case.  That  is 
hot 

H  -  L2(H 

<  f.g>  «  f  f(x)g(x)  d 

Jr  _ 

Then  <  Af,g>  “  <  f,Ag  >  .  This  is  real  symmetry  and  A  Vf  A*  i.e.,  A 

is  non-self  adjoint. 

Question.  When  is  the  Ansatz  correct?  Sufficient  condition:  AA*  -  A*A 
i.e.,  A  is  normal  .  Then  an  orthogpnal  basis  can  be  found  in  H  =  L^IO. 

Here  (1)  requires 

r  sinkt  |x-tj -|  t-yj  } 

ir  |£-t|  |t-*|  d°  t  =  0 

This  last  condition  can  be  shown  to  be  satisfied  by  the  cylinder  and  sphere 
but  not  for  the  ellipse  or  ellipsoid.  In  particular  then  any  EEM  theory  results 
which  use  the  Picard  process  and  are  used  to  construct  equivalent  circuits  arc 
suspect  in  general. 
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before  entering  into  what  is  known  in  the  case  when  the  operator  A  is 
not  r-.tTT3.il  the  relation  between  SEM  and  EEM  when  the  Picard  process  is  valid 
should  be  mentioned. 


In  [2.21  -Baum  discussed  the  matrix  case  and  showed  that  every  zero  of  Xfc) 
was  a  pole.  Kore  generally: 

The  Relation  between  SEM  and  EEM 

Theorem  (P.a.Tur.i.  The  poles  of  G(x,y,k)  are  zeros  of  the  eigenvalues 

l  eik| x-y| 

n '  "  *  J0 


(26) 

(27) 

where 


A-(k)  =  0  G-  =  4? 


(x,y,k)  **  G0(x-y)  -  £  GQ  (x,  s, k)U  (s,y , k) 
ik|x-y| 


d7 


3G_ 

<?n 


9  0 

Bo  *-  r TT— r— 1 —  p(y)do 

2n  2i*  I  x-y|  y 


U  ♦  E'J  =  2 


!fc 

3n 


If  the  operator  is  not  normal  the  situation  is  much  more  complicated  in 
general.  One  usually  has  to  contend  with  root  vectors  as  they  occur  in 
the  Jordan  normal  form.  The  simplest  example  of  their  occurrence  is  in  the 
matrix  solution  of  ordinary  differential  equations  with  repeated  roots.  The 
questions  of  when  are  the  root  vectors  complete,  when  do  they  form,  a  basis 
are  difficult  in  general .  There  is  one  case  when  there  is  a  simple  theorem 
concerning  co-pi  oneness,  namely  if  the  operator  is  dissipative.  An  operator 
A  is  said  to  ho  dissipative  if 
lr.  «•  A : ,  i  >  10. 

Many  of  t'.<.  operators  in  mathematical  physics  are  dissipative.  For 
example  the  free  space  Green's  function  is: 
r  _ ilti  x-y  | 


$(x)$(y)d3x  d3y 

kernel  implies  that 
0  . 

A  rigr.  :rou:.  preof  of  this  can  be  found  in  Dolph  (6.20).  * 

Fair.-,  (6.72)  has  established  a 'completeness  theorem  for  such  operators 
which  are  compart  and  nuclear. 

Before  stating  his  result  note  that  if  S  (A)  are  the  eigenvalues 

1/2  n  f? 

of  (A *A!  a  compact  operator  is  called  nuclear  if  )  S  (A)  <  “  . 

l'  N 

Theorem  1.  If  A  •-  P  +  N  where  P  is  positive  and  compact,  and  U  is 
dissipative  and  nuclear.  Then  the  toot  vectors  of  A  are  complete. 


j  <t|x-y| 
One  has  for  real  k 


<■  (y)dy 


Xr  < 


...  (!  sin  k | x-y 

I) 

iko  behav 

'  '  >  3  [l* 


and  the  delta  like  behavior  of  the 

Im  <  A), ; 


(x) I 2d3x  > 
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A  simple  pertinent  example  is  given  by 
silc|x-y| 


A* 


-£  f 


l*-y| 


$(y)da 


by  taking  P)  *=  -p-  f  and  =  (A-P)^. 

4’'  J  |x-y| 

More  information  on  root  vectors  and  basic  can  be  found  in  the 
reference  (6.34). 


Finally,  space  limitations  do  not  permit  detailed  discussion  of  many 
topics  important  to  the  further  development  of  this  subject.  These  include 
the  weak  perturbation  of  compact  operators  see  (6.42,  (6.55),  (6.75)  as 
well  as  variational  principles  (6.74)  and  papers  in  press  by  Ramm. 
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INTRODUCTION 

This  is  a  brief  summary  of  the  invited  talk  given  by  the 
author  at  the  Lexington  (November  1980)  meeting.  The  purpose  of 
this  paper  is  to  formulate  the  mathematical  problems  important 
for  the  SEM  and  EEM,  to  answer  several  basic  questions  and  to 
draw  attention  to  certain  unsolved  problems.  Some  r.ev;  results 
are  also  reported.  The  detailed  presentation  of  the  talk  was 
sent  to  the  Mathematical  Notes  (ed.  C.  E.  Baum)  and  submitted  fo; 
publication  in  the  J.  Math.  Anal.  Appl.  The  bibliography  is  not 
complete:  only  the  papers  in  which  the  results  mentioned  it.  thi, 

article  appeared  were  included  in  the  bibliography. 


1.  STATEMENT  OF  THE  EEM  AND  SEM 

Let  f  be  an  e-xterior  domain  with  a  smooth  closed  boundary  r 
D  be  the  corresponding  Interior  domain, 

G0  =  ~f  x'lf  i  ' »  r  =  M>  AS  =  Jp  G0(s,s' ,k)g(s')ds'  and 
u  =  Jj.  Gq(x,s  ,k)g(s)dx.  The  function  u  solves  the  problem 

(V2+k2)u  =  0  in  D,  u|p  =  f,  r(3u/3r-iku)  -*•  0  as  r  ♦  ~  (1) 

provided  that 


Ag  =  f  (2) 

If  one  uses  the  Laplace  transform  variable  p,  then  p  =  -ik,  and 
the  half  plan?  Rep  >  0  corresponds  to  the  half-plane  Imk  >  0. 
Engineers  (6.41]  -  (2.2)  tried  to  solve  (2)  by  the  formula 

g  —  Aj  c  j  f  j  ,  where,  A  f  j  —  I  j  1  j  >  I  ^  ^  I  ^  I  Ap  |  —  *  *  *  9nd 

f  =  This  can  be  done  if  A  =  A*  Is  selTadjoint  on 

H  -  L2(f).  The  operator  A  In(2)  is  nonsel fadjolnt .  Therefore: 
1)  it  may  have  r.o  eigenvectors  (e.g.  Ag  =  Jggdx  on  H  =  L2[0,1]) 
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2)  it  nay  have  not  dnly  eigenvectors  but  also  root  vectors  16.69), 

fc.63, 

3)  it  is -an  open  question  whether  one  can  expand  an  arbitrary 
function  fcH  in  the  series  of  eigenvectors  and  root  vectors  of  A. 
Of  course  one  is  interested  in  the  rate  of  convergence  of  the 
series  in  eigen  ar.d  root  vectors  and  in  algorithms  for  calcula¬ 
tion  of  the  roc:  vectors  and  eigenvalues  of  A.  The  outlined 
method  ( EE>: )  r.as  the  following  merits:  1)  instead  of  problem 

(1)  with  a  continuous  spectrum  in  the  unbounded  domain  we  con¬ 
sider  problem.  (2)  with  a  discrete  spectrum  on  the  compact  mani¬ 
fold  r,  2)  the  resonance  properties  can  be  conveniently  studied 
by  the  ESK.  A  mathematical  study  of  the  EEM  was  originated  in 
[6.72),  [6.73],  [5. 71]. 

In  order  to  describe  SEM  consider  the  problem 

utt  =  v?u  in  n»  u! r  =  °»  uU=o  =  °*  ut|ts£)  “  f  (3) 

The  solution  of  this  problem  takes  the  form 

u  *  (2it)"*  exp(-ikt)v(x,k)dk,  v  *=  JfiG(x,y,k)f  dy  (*l) 

where  G  is  the  Green  function  for  problem  (1), 

G  *  G0  -  fr  G0(x,s,k)uds,  p  *  — - 

s 

3Gn  3Gn 

[ I  +  T ( k )  ] u  =  ?  jf;  T(k)p  =  2/pgor  V  as’ 

Wo  assume  that  fcC^Cil).  From  (5),  (6)  it  follows  that  G  is 

ir.ite-neromorphic  in  k.  This  means  that  G  is  mesomorphic  on 
he  whole  complex  plane  k  and  its  Laurent  coefficients  are  de- 
e.nerate  kernels  (finite  rank  operators  on  H) .  If  ft  C  then 
ic  analytic  in  Im  k  >  0.  Thus  v  is  meror.orphic  in  k  and  analy¬ 
se  in  Im  k  _>  0 .  Let  us  assume  that 

!  v !  <  c(b:  (l  +  |k|  )~3,  a  >  1/2,  J Re  k|  -  Im  k  =  b  (7) 

where  b  is  ar.  arbitrary  const; 

I  Im  X ,  j  ♦  ■  as  J  +  »,  |  Im  kj  |  <  |  Im  V.?\  <  . . .  (8) 

where  k,  are  the  poles  of  v. 

Note  that  { '7 )  =£>  ( S ) .  From  (7),  (8)  it  follows  (by  moving 
the  contour  of  integration  In  (U)  down)  that 

u(x,t)  -  c  i  ( x  ,t)exp(-ikjt)  +  o(exp(-|lm  kH|t)),  t  •»  +  «°  (9) 

Here  c.(>.,t)  expC-ik.t)  =  F.es  v(x,k)  exp'(-ik.t)  at  k  =  k.; 

TV;  .  —  1  **  « 

c.(y.,t)  =•  C  ( t  -  ),  where  m,  is  the  order  of  the  pole  fc..  Thus 

J  J  J 


(5) 

(6) 
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we  see  that  (7),  (8)  and  the  meromorphic  character  of  v  are  suf¬ 
ficient  for  the  SEM  of  the  form  (9)  (asymptotic  SEM ) .  It  is  an 
open  question  if  the  series 

u(x,t.)  =  Cj(x.t)  expC-lkjt)  (10) 

converges.  The  validity  of  the  EEM  was  discussed  in  (6.691,  (6.681. 


2.  COMPLEX  POLES  OF  CREEK'S  FUNCTIONS 

We  saw  ir.  Section  1  that  complex  poles  kj  are  important.  It 

is  interesting  to  answer  the  following  questions:  1)  how  does 
one  calculate  the  poles?  2)  are  the  poles  simple?  3)  do  the 
poles  depend  continuously  on  the  scatterer?  4)  Can  one  identify 
the  scatterer  from  the  knowledge  of  complex  poles?  5)  what  can 
be  said  about  location  of  the  poles  and  asymptotic  behavior  of 
the  large  poles  nearest  to  the  real  axis?  6)  are  there  any  mono- 
tor.icity  cr  other  features  in  the  behavior  of  the  purely  imaginary- 
poles?  7)  What  are  the  properties  of  the  resonant  states  (natural 
modes  corresponding  to  the  complex  poles)?  8)  What  is  the  rela¬ 
tionship  between  the  poles  and  the  eigenvalues  used  in  the  EEM? 

We  give  some  answers  to  the  above  questions.  Three  different 
methods  for  calculation  of  the  complex  poles  were  giver,  in  [6.71-2], 

[6.68]  and  [6.7-1] .  The  first  method  is  most  general.  It  reduces  the 
problem  to  calculation  of  the  values  k^  at  which  a  certain  opera¬ 
tor  of  the  type  I  +  T(k),  where  T(k)  is  a  compact  analytic  opera¬ 
tor  function,  is  not  invertible.  These  kj  can  be  found  by  a 

projection  method.  The  method  is  described  in  [6.71-2]  (see  also 

[6.68] ).  Its  convergence  is  proved  16.71]  The  second  method  is  a 

2 

variational  principle  for  complex  poles:  kj  are  the  stationary- 
values  of  the  functional 

K(u)  =  <Vu,  Vu>/<u,u»,  where  <u,v>  =  lira  J  exp ( -e | x | ln| x | )uvdx 

e-+0  (10’) 

and  the  integral  is  taken  over  ft.  In  [6.74]  a  certain  system  of 
test  functions  was  suggested  but  the  rigorous  justification  of 
the  numerical  procedure  given  in  [6.74]  is  an  open  problem.  In  [6.68] 
a  variational  principle  for  the  spectrum  of  compact  nonselfadjoint 
operators  was  given.  In  [6.71] it  was  proved  that  the  complex  poles 
of  the  Green's  functions  are  the  complex  zeros  of  the  eigenvalues 
of  certain  integral  operators.  This  gives  the  third  nethod  of 
calculation  of  the  poles:  first, one  calculates  the  eigenvalues, 
then  one  looks  for  their  zeros.  No  numerical  results  are  known 
for  the  third  method.  It  would  be  interesting  to  make  numerical 
experiments  and  to  compare  all  tl)e  three  methods. 

It  is  an  open  question  whether  the  poles  are  simple.  In  [6.71] 
it  was  proved  that  the  poles  are  simple  if  the  surface  is  of  such 
shape  that  the  operator  A  in  (2)  is  normal,  i.e.  AA*  =  A*A.  In 
16.73]  It  was  proved  that  this  is  so  if  f  is  a  sphere  or  a  straight 
line  (linear  antenna).  Recently  the  author  gave  a  simple  example 
of  a  multiple  pole  in  the  problem  with  third  boundary  condition: 
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if  (V?+k2)  u  «  0  ir.  r  *  |x|  _>  1,  3u/3r  -  2u  =  cos  8  on  |x|  *>  1, 

r(-|p  -iku)  -»  0  as  r  ■*  •»,  then  k  =  -21  is  a  pole  of  order  2  of 

u(x,k).  Gersrieally  multiple  poles  are  exceptions  because  small 
perti: rbat io' s  of  the  shape  of  the  scatterer  can  destroy  multiple 
poles.  On  t h ;  other  hand,  since  the  poles  depend  continuously  on 
f  (see  [6.63’  for  precise  definitions  and  proofs)  it  seems  possible 
that  by  continuous  variation  of  T  one  can  make  a  multiple  pole 
out  of  2  sirrle  poles  by  merging.  Nevertheless,  no  proof  is  known 
that  the  Green’s  function  of  the  exterior  Dirichlet  Laplacian  has 
multiple  pole.-  for  some  f. 


Vie  have  already  mentioned  that  the  poles  depend  continuously 
on  T.  It  is  r.ot  known  whether  the  set  of  complex  poles  deter¬ 
mines  the  scatterer  uniquely.  A  discussion  of  this  question  is 
in  [6.76)ar.d  [  *].  Some  information  on  location  of  the  poles  is 
available:  in  [6.70] it  was  proved  that  the  domain  (Im  k  <  0, 

| Im  k|  <  alog|Re  k|  +  b,  a  >  0}  is  free  from  the  complex  poles 
of  the  Green’s  function  of  the  Schrodinger  operator  with  a  com¬ 
pactly  supported  potentie1 ;  in  [6.54]  a  similar  result  was  proved 
for  the  poles  of  the  Green's  function  of  the  exterior  Dirichlet 
Laplacian;  in  [6.5)  some  heuristic  arguments  are  given  to  show 

that  the  domain  {In  k  <  C,  |  Im  k|  <  a|Re  k|*^  +  b,  a  >  0}  is 
free  from  the  poles  of  the  Green's  function  of  the  exterior 
Dirichlet  and  Neumann  Laplacians  provided  that  r  is  strictly  con¬ 
vex  ar.d  smooth;  if  f  is  not  smooth  (say,  r  is  a  polygon)  then  there 
exists  a  series  cf  poles  kj  such  that 

| In  V: j  j  =  0(ics|Re  k,|)  as  J  ■*  <»  (6.6). 

In  [6.53]  i;  was  proved  that  there  exist  infinitely  many  purely 
imaginary  poles  of  the  Green's  functions  of  the  exterior  Dirichlet 

or  Neumann  I.ar'.  acian  a.ni 

cR^  <  1  im  inf  y~2  H(y)  <  lim  sup  y~2  N(y)  <  cRp 

V  -*  -*•  y->w 


where  c  *=  1.135370 . N(y)  is  the  number  of  purely  imaginary 

poles  with  | In  k.|  <  y,  the  obstacle  is  star-shaped  (this  means 
J 

that,  all  points  c  f  T  can  be  seen  from  a  point  in  D)  and  Rj ,  R^ 

are  the  radii  o”  spheres  inscribed  in  and  circumscribing  D, 
respectively.  It  is  pointed  out  in  [  *]  that  if  D-,  =  qD,  ,  q  >  1 

M  \  (11.  (?)  £  A 

then  yj  =  qy]“',  where  -iyj  '  (-iyj  >  are  the  poles  of  the 
Green's  function  cf  the  exterior  Dirichlet  Laplacian  in 

ilj  •=  R  -\ D ^  ,  i  =  1,  2,  where  R^\D  denotes  the  complement 

to  D  in  R^.  Therefore  in  this  case  N-,(y)  >  N.(y)  and 
(1 )  (?)  ( n  £  ~  1 

y^  >  y. '  ,  where  are  thd  moduli  of  the  purely  imaginary 

poles  with  minimal  moduli.  In  [6.53)  Theorem  3-5  on  p.  751  says 
that  !•!  (y)  <  ‘.^(y).  This  statement  contradicts:  1)  the:  above 

argument. ,  ar.d  2)  the  case  when  D^  and  D^  are  concentric  balls  and 

one  c r.  calculate  .’.'j(y)  and  N^(y)  for  y  >>  1  and  verify  that 
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>  I.^Cy).  The  argument  inl6.53)  can  be  used  if  the  assumption 
0  C  0^  is  replaced  by  the  assumption  0  3  0s>  We  mention  this 

because  in  the  literature  one  can  find  references  and  citations 
of  Theorem  3-5  from  [6.53]  in  its  wrong  form.  Using  arguments  from 
[6.53]  and  assuming  that  Dj  ,  j  =  l,  2  are  star-shaped  and  that 

D.^  C  D-,  C  D^,  one  can  see  that  N^y)  £  N^ty)  <_  N^(y).  here  we 

used  the  corrected  version  of  Theorem  3-5  from  [6.53]:  if  D,  C  D„ 

1  e 

and  is  star-shaped,  then  N ^ { y )  <_  N^(y).  This  theorem  is  actu¬ 
ally  proved  in  [6.53]  so  that  the  misstatement  of  Theorem  3-5  in 
[6.53]  is  just  a  misprint. 

Concerning  the  behavior  of  the  resonant  states,  that  is  the 
solutions  of  the  homogeneous  problem  (1)  for  k  =  a  -  iy,  y  >  0, 
f (x)  =  0),  satisfying  the  asymptotic  condition 

u  =  I x | -1  exp(ik|x| >lj_0  lxi~Jfj»fj  c  fj  v  n , y ) ,  n  =  x  ♦  |x|_1,  (11) 

at  infinity,  one  can  prove  the  following  proposition:  if 
u  exp(-y  |  x  |  ) !  x  j  -»  0  as  |x|  -»  «•  then  u  =  0.  From  this  it  follows 
that  the  resonant  states  (scattering  modes)  corresponding  to  a 
complex  pole  k  =  a  -  iy  grow  at  infinity  exactly  as 

0(exp (y | x |  )  | x J _A .  See  also  [6.43]  Theorem  3.  The  relationship 
between  SEX  and  EEM  is  given  in  the  following  proposition ([6.71-2] , 
[6.68]):  the  se*  of  the  complex  poles  of  the  Green's  function  of 
the  exterior  Dr'richlet  Laplacian  coincide  with  the  set  of  complex 
zeros  of  the  eigenvalues  *nOO  of  the  operator  A  defined  in  (2). 

It  is  not  known  at  this  time  whether  the  order  of  a  pole  can  be 
calculated  from  the  multiplicity  of  zeros.  One  can  construct 
other  operators  with  the  eigenvalues  vanishing  at  the  complex 
poles  [6.68] . 


3.  "ORTHOGONALITY"  OF  THE  EIGENMODES  AMD  RESONANT  STATES 


By  eigennoies  (EM)  we  mean  the  eigenfunctions  of  the  opera¬ 
tor  A  defined  in  (2).  This  is  a  nonselfadjoint  operator  on 

H  =  L?(D  with  the  property  [Af,g]  =  Lf,Ag],  where  [f,g]  =  (f,g)  = 
Jj.fg  qs  ,  (•,*)  is  the  inner  product  in  L  (D,  the  bar  denotes 
complex  conjugation.  Suppose  that  Afj  =  Ajfj,  [fj.fjl  K  0,  j  - 
1,  2,...  arid  the  set  { f  j }  forms  a  basis  of  H.  Then  any  fcH  can 
be  represented  as  f 


r® 

ij  =  l 


Cj  fj  and  Cj  =  l  f ,  fj  ] . 


This  can  be 


proved  exactly  as  in  the  case  of  orthogonal  Fourier  series  if  one 
takes  into  account  that  [fj.f^]  =  0  for  j  ^  m.  The  last  formula 

follows  from  the  identity  0  =  [Afj,fm]  -•  [fj»Afml  =  ( A j  — ) 

If  A.  =  A  one  can  choose  f,,  f  so  that 
J  m  J  m 

[f,,f  ]  =  0  for  j  ^  m.  Thus  the  coefficients  in  the  EEM  can  be 

J  *■! 

easily  calculated.  If  the  root  vectors  are  present  the  formulas 


[fj.rj  ir  x  x„. 
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for  the  coefficients  in  root  vectors  can  also  be  calculated  ex¬ 
plicitly  [  *  ]. 

"Orthogonality"  of  the  resonant  states  corresponding  to  dif¬ 
ferent  complex  poles  ,  k2  holds  in  the  following  sense: 

<u(x,k^),  u(x,k2)>  =  0,  where  <•,•>  is  defined  in  (10‘)  (See  {6.74] 

and  [*  ]  for  detailed). 

4.  /JO: 'SMOOTH  BOUNDARIES 

The  usual  proof  of  the  meromorphic  nature  of  the  Green's 
function  of  the  exterior  Laplacian  requires  smoothness  of  T. 
Indeed, it  is  based  on  the  integral  equation  (6)  and  on  the  theorem 

about  the  neronorphic  nature  of  the  operator  (I+T(k))_'*‘  [6.80-3].  if 
this  theorem  it  is  assumed  that  T(k)  is  a  compact  operator  func¬ 
tion  analytic  in  k.  If  the  surface  T  has  edges  or  conical  points, 
the  operator  T(k)  in  (6)  is  no  longer  compact.  Nevertheless  the 
theory  is  still  valid  provided  that  there  are  no  cusps  on  r. 

This  follows  from  the  proposition  (see  [  *  ]  for  details):  if 
T(k)  =  T  +  Q(k)  is  an  operator  function  on  a  Hilbert  space  H, 
where  Q(k)  is  analytic  in  k  for  keA,  where  A  is  a  connected  open 
set  in  the  complex  plane,  |T|ess  <  1  and  I  +  T(k)  is  invertible 

at  some  point,  then  (I+T(k))~1  is  finite  meromorphic  in  A,  (finite 
merc.-.orphie  means  that  the  Laurent  coefficients  are  operators  of 
finite  rank),  by  j T | ^ s._  we  mean  inf  |ff-K||,  where  K  runs  through 

the  set  of  all  compact  operators  on  H. 

It  is  known  [6.12]  that  J  T  (0 )  { esti  <  1  provided  that,  there  are' 

no  cusps  or:  f.  One  can  now  apply,  the  above  proposition  and  con¬ 
clude  that  u  i:.  (5)  (ar.d  therefore  G;  see  (0))  is  meromorphic  and 
its  Laurent  coefficients  are  degenerate  kernels. 


5.  EXAMFLSS,  COMMENTS 

1.  A  symmetric  (with  respect  to  the  form  [f,g]  defined  in  section 
3)  r.onselfadj  oint  operator  can  have  root  vectors.  Example: 

1  ■*  •*  *“1 
A  =  (i_t),  [x,Vj  =  *1^1  +  x2y2'  ( A-XI )  has  a  pole  of  order  2 

at  ).  =  0.  The  corresponding  eigenvector  is  (!)  and  the  root 

vector  is  (1^*). 

2.  The  fact  that  the  algebraic  problen  to  which  an  original  inte¬ 
gral  equation  was  reduced  (e.g.  by  a  projection  method)  has  eigen¬ 
values  does  net  guarantee  that  the  original  equation  has.  (See 
[6.68],  [6. 69]  and]*  ]  for  details  and  sufficient  conditions  under 
which  the  eigenvalues  of  the  algebraic  problem  converge  to  the 
eigenvalues  of  the  original  problem.) 

__  —  1 

3.  The  operator  (1+T(k))  '  cat;  have  multiple  poles  and  be  dia¬ 
gonal  i cable  (i.e.  T(k)  has  no  root  vectors). 
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ixar.pl e :  ?Ck)  =  (~l  +  k?  °-V  (I+T(k))-1  =  /k  ~  0 

\  o  >c  y  \o  d+k2)-1/ 


Exact 

k  =  0  is  a  r :.le  of  order  2. 


There  exists  an  operator  with  the  root  system  which  forms  a 
basis  of  h'  but  under  a  different  choice  of  the  root  vectors  the 
root  system  of  this  operator  does  not  form  a  basis  of  H  (see  [  *  ] 
for  an  example). 


t).  The  set  of  complex  poles  of  the  Green's  function  of  the 
Schrodir.ger  operator  does  not  define  uniquely  the  potential  if 
there  are  bound  states  [  *  ]. 

6.  If  z  is  a  complex  pole  of  order  m  of  (I+T(k))~\  T(k,e)  is 
compact  and  analytic  in  k  and  e  for  {|k-z|  _<  a,  |e|  <  b}  and 

T(k,0)  =  T(k)  then  the  poles  z(s)  of  ( I+T(k ,e ) ) can  have  a 
branch  point  at  e  =  0  and  ord  z(e)  £  m.  Moreover  z(e)  can  be 

1/r 

represented  by  Puiseux  series,  i.e.  by  a  series  ir.  powers  of  c  , 
where  r  Is  some  integer  (see  [*  ]  for  details). 


7.  The  multiplicity  of  the  complex  poles  is  not  equal  to  the 
order  of  zeros  of  eigenvalues,  generally  speaking. 

It  was  proved  ir,  (6.72](see  also(6.68J)  that  the  set  of  complex 
poles  coincide  with  the  set  of  complex  zeros  of  the  eigenvalues 
of  certain  integral  equations.  In  the  case  we  are  concerned  with 
in  this  paper  cr.e  can  have  in  mind  the  eigenvalues  of  the  equa¬ 
tion  [1+T(k)]  u.  =  Xj(k)Uj,  J  =  1,  2,....  It  was  an  open  question 

whether  the  orders  of  the  zeros  of  Xj(k)  are  equal  to  the  multi- 

plicities  of  the  corresponding  poles.  We  show  by  presenting  an 
example  that  this  is  not  so  In  general.  Let  us  take  as  I  +  T(k) 
a  finite  dimensional  operator  with  the  following  matrix 

A(k)  =  />(k)  1  We  have  X  (k)  «  X(k),  A-1(k)  =  (X  1(k)  ~X~!(k)V 

\  0  X(k )/  J  \  0  X"1 (k)/ 

If  X(z)  -  0  and  m  is  the  order  of  the  zero,  then  z  is  the  pole 

of  A~‘(k)  of  multiplicity  2m.  It  is  clear  from  this  example  that 
the  order  of  zeros  of  the  eigenvalues  will  coincide  with  the 
multiplicity  of  the  correspond ing  poles  iff  A(k)  is  diagonalizable , 
that  is  A(k)  has  no  root  vectors.  This  example  is  sufficiently 
general  because  for  a  compact  T  the  eigenvalues  \j  t  -1  have 

finite  algebraic  multlpl icites  and  the  corresponding  root  spaces 
reduce  I  +  T(k),  so  that  in  the  root  spaces  I  +  T(k)  is  a  matrix 

opera tor. 

8.  U-'ing  the  ideas  given  in  (6.68)  the  author  proved  convergence 
of  th --matrix  approach  in  scattering  theory,  widely  used  in 

p  ra  c  t  *  c  . 
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9.  A  variational  principle  for  complex  poles 


In  section  3  it  was  mentioned  that  the  complex  poles  of  the 
Green  function  occur  at  the  complex  points  k  at  which  the  homo¬ 
geneous  equation  (2)  has  a  nontrivial  solution.  Let  H  denote 

the  Sobolev  space  V.'^(r),  and  |f|^  denote  the  norm  in  H^.  Consider 
the  variational  principle  F(f)  =  |Af|^  =  min,  | f | q  -  1.  If  {fj} 


/  «  n 

is  a  basis  of  H  =  H0,  and  f'  '  =  £  c.f,,  then  the  problem 

J=1  J  J 


F( 


(f(n))  =  mir.,  |f^n;|n  =  1  yields:  l  a 

n  m=l  J 


.(n) 


0Ocr 


(n) 


°,  l<  3  in, 

be  the  complex 


l  | c  .  |  >  0.  Thus  (*)  det  a  (k)  =  0.  Let  k 

j_l  J  jm  _ 

roots  of  ( * ) .  Then  it  can  be  proved  that  the  set  of  the  complex 

limit  points  {k  )  of  the  set  { k ^ n ^ )  coincides  with  the  set  of  the 
s  s 

complex  poles  of  the  Green  function.  This  a  new  result.  The 
functional  F(f)  is  real  valued  in  contrast  with  the  functional 
K(u)  in  (10'). 


Problems 

1)  Is  it  true  that  the  root  systems  of  A(k),  T(k)  form  a  Riesz 
basis  of  H?  It  is  proved  that  these  systems  form  a  Riesz 
basis  with  brackets  (see  [6.68] for  a  proof  and  definitions). 
The  author  thinks  that  the  answer  is  no. 

2)  Is  there  a  relation  between  the  order  of  a  complex,  pole  and 
the  multiplicity  of  the  zeros  of  A. n ( k )  ? 

3)  Can  the  scatterer  be  uniquely  identified  by  the  set  of  com¬ 
plex.  poles  of  the  corresponding  Green's  function? 

*0  Prove  that  there  are  infinitely  many  complex  poles  kj  with 
Re  k,  *  0  (in  diffraction  problems  and  noncentral  potential 

v 

scattering) . 

5)  Are  the  complex  poles  of  the  Green's  function  of  the  exterior 
Dlrichle;  or  Neumann  Laplacian  simple? 

6)  Make  numerical  experiments  in  the  calculation  of  the  complex 
pcies . 

7)  Prove  convergence  of  the  numerical  procedure  for  calculation 
of  the  complex  polos  suggested  in  [6.74]. 

8)  Find  a  theoretical  approach  optimal  in  some  sense  to  approxi¬ 
mate  a  function  f(t)  by  the  functions  of  the  form 

**  ^  , 

N  J  IH-*  1 

1 ^  exp(-ikj t ) t  ~  cmj .  Here  the  numbers  m  ^  ,  kj 

are  to  be  found  so  that  f^  will  approximate  f(t)  in  some 
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optimal  way.  Currently  some  methods  (e.g.  Prony  method)  are 
used  It;  practice,  but  they  are  not  optimal.  This  problem 
seems  to  be  of  general  interest  (optimal  harmonic  analysis 
in  complex  domain). 

9)  When  can  SEM  in  the  form  of  (10)  be  justified? 


Conclusion 

We  hope  that  it  was  shown  in  this  paper  that: 

1)  EEM  is  justified  (in  the  generalized  form  of  expansion  in 
root  vectors). 

2)  SEM  is  justified  in  the  asymptotic  form  (9). 

3)  Numerical  projection  method  for  calculation  of  the  complex 
poles  is  Justified. 

4)  There  are  many  interesting  and  difficult  open  problems  in 
the  field. 

5)  Numerical  results  and  experiments  are  desirable. 

Reference 

(  *)  Ramm,  A.G.  ,  Mathematical  foundations  of  the  singularity 

and  eigenmode  expansion  methods,  J.  Math.  Anal.  Appl.  (1981). 
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1.  INTRODUCTION 

At  the  international  symposium  on  wave  scattering’ 
most  of  the  speakers  pointed  out  that  the  7-matrix  scheme 
needed  a  justification,  i:s  convergence  was  not  proved.  In 
this  paper  a  proof  of  car-  ergence  is  given.  This  proof  also 
clarifies  another  basic  .ju.:st:or!,  namely,  convergence  of  the 
variational  method  of  finding  stationary  points  of  function¬ 
als.6  Many  physical  problems  are  formulated  as  the  prob¬ 
lems  of  finding  stationary  points  and/or  stationary  values  of 
some  functionals,  and  these  points  are  not  extremal.  A  neces¬ 
sary  and  sufficient  condition  for  a  stationary  principle  to  be 
extremal  is  given  in  Ref.  2,  p.  90.  The  standard  7-matrix 
approach  is  described  in  Ref.  1,  pp  64.  The  principal  differ¬ 
ence  between  the  standard  and  our  approach  is  as  follows.  In 
the  standard  approach  the  scattered  field  is  represented  as 
the  series  in  the  outgoing  spherical  waves  and  the  coeffi¬ 
cients  of  the  series  are  found  from  a  linear  system.  One  as¬ 
sumes  that  the  series  converges  on  7  (the  Rayleigh  hypoth¬ 
esis)  which  is  not  true  in  general.  In  cur  approach  one  uses  a 
basis  of  L  3(7)  and  no  difficulties  with  convergence  arise. 

Let  us  describe  a  modified  7-matrix  approach  to  the 
problem 

(_VJ -*’-):/  =  /  in  fl,  A>  0  (1) 

u|r=0,  r\3u/dr  —  iku, — 0  as  r  =  |xj— ►  <»  (2) 

Here/}  isan  exterior  domain  with  a  smooth  dosed  boundary 
r and  D  —  R}\/2  is  a  bounded  domain.  From  the  Green 
formula  it  follows  that 


u{x)  =  c{x)  =  J  g{  xj)h  (s)ds,  xeP. 

(3) 

0  =  ,(*)  —  f  si  xj]h  xsD 

(4) 

•W-  (t/o. 

Jn  4rr|je  — 

A -*L. 
dN 

(5) 

JVis  the  unit  normal  to  7  directed  into/?.  If  A  is  found,  then 
u[x)  can  be  found  from  (3).  Let  us  rewrite  Eq.  (4)  as 


Ah  =  f(r).  Ah  ~  J  jfaj’)hds‘,  ssT.  (6) 

Let  [<f>j \,j  =  1,2,. . .  be  a  basis  of  where 
H  —  H0—L\r),  H,  —  W 5(7)  are  the  Sobolev  spaces,14 

K  -  2  <A>  W 

j- i 

i  V-=V  (8) 

m  m  I 


' ■Supported  by  AFOSR  8TC2CW. 


where 

ajm=[A6„A).  Vj  —  [u.tfij),  (fv)  =  (/,«,)„.  (9) 

Let  (f^)  be  the  inverse  matrix  a\~  ",  l<y,m<rt.  Then  cJt 
1  <  j<n,  can  be  calculated  if  Vj,  1  <  j  <n,  arc  given.  From  for¬ 
mula  (7)  one  calculates  A. ,  and  from  formula  (3)  with  A  =  A. 
one  calculates  the  approximate  solution  u„  to  problem  (!}- 
(2).  The  problem  is  to  prove  that  (i)  for  sufficiently  large  n  the 
matrix  ctj„  in  (8)  is  invertible,  (ii)  ||A„  —  h  J( — »0  as  n— »oo, 
where  |j-||  is  the  norm  in  H.  Actually,  convergence  will  be 
proved  in  ,  where  q  depends  on  /and  on  the  smoothness  of 

7.  Let  us  assume  for  simplicity  that  7C  C  ",  Then  q  depends 
on  the  smoothness  of/if  meas(7n  supp f)  >  0  (supp /is  the 
support  of /)  and  q  is  arbitrary  if  dist(supp/,7)>0. 

The  basic  idea  of  the  proof  is  very  simple  and  is  given  in 
Sec.  2.  In  Sec.  3  some  technical  details  are  provided. 


2.  MAIN  RESULT 

Theorem  1:  System  (S)  is  uniquely  solvable  for  suffi¬ 
ciently  large  n  and  ||A„  -  h  ||-^)  as  n—oo  [without  loss  of 
generality  we  assume  that  the  operator/  +  7  (A)  is  invertible; 
see  the  proof  below,  n.I  of  Sec.  3]. 

Proof-.  The  basic  idea  is  to  factorize  A  in  (6)  as 
A  =  At>[  /  +  7(A )),  where  A(>  =  A  (0).  T{k )  =  A  0~  1 
X[A  [k)~  A  (0)].TheoperatoM„>0isabijectionof//ff  onto 
Hn  +  , ,  while  T  (A )  is  compact  in  //,  for  any  q  (see,  e.g.,  Ref.  2. 
p.  287).  The  system  (8)  can  be  written  as 
[At,l  /  +  7}A„,f$;)  =  (w,^),  1  <,j<,n.  Since  ^o>0  the  form 
[Au  u,f )  is  a  scalar  product  which  we  denote  by 
[u,f]  =  [A0  u.f).  This  scalar  product  generates  a  norm 
which  is  equivalent  to  the  norm  in 
H_  ,/j.  This  follows  from  the  fact  that  A0  is  a  pseudodifTeren- 
tial  elliptic  operator  of  order  —  1  and  therefore 
ord  A  =  —  J.  Thus  (8)  is  of  the  form 
[(  /  +  T)h„,6j  ]  =  (v,dt),\<,j<n.  Lttw~A  0~  'u.  Then 
=  [w.<f>j  Jand  A„  -f  7„  Thm  =  P,w,  where  7„  is  the 
orthoprojection  in  H_ ,/:  on  the  linear  span  of  .  .,d„  |- 
The  operator  /  +  T (A ),  A  >  0,  can  be  assumed  invertible  (this 
will  be  shown  in  Sec.  3),  and  7 (A )  is  compact  on  //9  with 
arbitrary  q,  —  »  <q  <  m .  Therefore  Jj(f  —  PJTH,— *0  as 
n—*  oo ,  and  the  norm  is  the  norm  of  operators  on  //f  (this  will 
be  explained  in  Sec.  3).  Thus  I  -r  P.  T  =  I  +  T  —  (/  —  P„)T 
is  invertible  for  sufficiently  large  n.  This  means  that  system 
(8)  is  uniquely  solvable  for  sufficiently  large  n.  Furthermore, 

h-hm=\l  +T)~'w-V  +  PmT)P,w 

_  =  B  \I  —  {I  —  P[m'TB)~'P ,  ]tu,  (JO) 
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where  B  =  [f  +  T)~',  Plm>  =  J  —  Pm.  Thus 

<C|  110-11,1!^.  “••II,.  C.C|  =const>0.  (11) 

where  Q, .  =  P'"TB.  || Qm  ||f-0  as  If  ||«0, 
-IM.-'HI,  <  co  then  (11)  shows  that  ||A,  —  h  ||f— *0  as 
n— »oo  and  the  rate  of  convergence  is  given  by  the  rate  of 
decay  of  the  magnitude  ||.P'',|7'||,  asn— »oo.  In  order  that 
\\A  o’  <  oc  it  is  necessary  and  sufficient  that  ee//,  „  1/2 . 
This  is  so  if  /£//,  _  ,  (-0  j  because  in  this  case  +  ,  (/? )  and 
its  trace  u|re//f  „  .n.  Our  argument  shows  that  if/eZ.  2(/2 ) 
=  H[fl)  the  smoothness  of  !/|r  is  even  higher  than  we  need. 
This  completes  the  proof.  Theory  of  the  H f  spaces  and  the 
trace  theorems  can  be  found  in  Rtf.  4. 


3.  ADDITIONAL  DETAILS 

(1)  Let  us  show  first  that  /  +  T(k ),  k  >  0,  is  invertible. 
Since  T is  compact,  it  is  sufficient  to  show  that  the  nullity  of 
this  operator  is  trivial.  If  (/  —  T[k  )]/>  =  0  then  A  ( k  )h 

—  Aq[I+  T[k)\h  =  0.  Therefore  the  function 
u(x)  =  /r  g(.x,  s)kds  solves  the  homogeneous  problem  (1>- 
(2).  It  is  well  known  that  the  solution  of  { 1)— (2)  is  unique. 
Thus  u(.x)==0  .  n  FI.  If  k  5  is  not  an  eigenvalue  of  the  Dirichlet 
Laplacian  in  D  then  y(x)==0  in  D,  and  from  thejump  relation 
for  the  normal  derivative  of  u  one  derives  that  h  =  0.  If  k  2  is 
the  eigenvalue  of  the  Dirichlet  Lap'acian  in  Z>,  then  the  argu¬ 
ment  is  the  same  but  instead  ofgi.t,  y,  k )  in  (3)-(5)  one  should 
use  the  Green  function  g,\X,y,  k )  cf  the  exterior  of  a  small 
ball  B,  C D.  This  ball  is  so  chosen  that  k  -  is  not  an  eigenvalue 
of  the  Dirichlet  Laplacian  in  D,  =  D  \B,.  Obviously  such  a 
ball  can  be  found  (there  are  infinitely  many  such  balls). 

Remark  1:  The  idea  of  applying  g, (at,  y,  k  )  in  order  to 
deal  with  the  case  when  k  '  is  an  eigenvalue  of  the  interior 
problem  was  used  in  Ref.  3. 

(2)  Let  us  show  that  ||  P’"T  (j—  0  asn—  «.  Since  [6, )  is 
a  basis,  one  has  |[  P,’'f\\-*0  as  r.~*x  for  any  /e/Z.  Since  7" is 
compact  it  can  be  written  as  T  —  Ts  -j-  dN,  where 

II  ||  <fy.  f  >  — 0  as  <V—  o: .  and  Ty  is  finite  dimensional: 

TN  f  —  Ijv_  ,  (/,  v,  .  Clearly  it  is  sufficient  to  prove  that 
||  />,"TAf/||<<5„  ||/ ji,  where  ,5„--»0  as  n—x.  One  has 
||  P'"Tyf\\<  2;v„ ,  I'  P'"*jlf,  d-)|! 

<11/1!  Z;.,  ||  $t\\  |i  <5„  ||/||,  where  5„— *0  as  n~*x 

because  ||  Z>,'"u;/|j-— 0  as  /j— *  x,  l<y<;V.  In  this  argument 
||  •  ||  can  denote  any  norm.  What  is  essential  is  that  P{n'— »0 
Strongly.  In  particular  cr.ecan  use  the  norm  of//_,/2  pro¬ 
vided  that  the  system  [  <3,  |  forms  a  basis  of  Note  that 
//_,/jO//,s°that  if  the  system  \  d>,\  forms  a  basis  of 
then  every  element/e1/  can  be  represented  in  the  form 
/  =  2“_  ,  Cjdj,  where  the  series  converges  in  H_U2 ■  It  does 
not  converge  in  //.  generally  speaking,  but  there  exist  bases 
such  that  if/etf  then  the  above  series  converges  in  H.  For 
example,  such  a  basis  is  the  basis  consisting  of  the  eigenfunc¬ 
tions  of  the  operator  <40  (see  also  Lemma  1  below).  If  [6, )  is  a 
basisof//,  then  |  A  '06J  j  is  a  basis  of //f  +  ..This  follows  from 
the  fact  that  A  '0  is  a  bijection  of  Hf  onto  ZZf  +  ,[A  '0  is  an 
elliptic  pseudodifierentia!  operator  of  order  —  s).  Since  |/|f 
<1/1,  for  q  <s,  it  is  clear  that  if  the  series  Z/_  ,  cy<Sy  con¬ 


verges  to /in  H  =  H0  it  converges  to /in  //_  ,/2.  It  is  conve¬ 
nient  to  have  a  system  {<Jyj  which  forms  a  basis  in  any  of  ZZf 
and  if /eZZ,  the  series  f—  2/_  i  c,d2  converges  in  .  For 
example,  if  S 1  is  the  unit  circle  and  //,  =  ZZf  (S ')  then  the 
system  |  exp{inx)/\/2rr |  forms  a  basis  of  ZZf  for  any  q.  The 
same  property  has  the  system  |  |  of  the  eigenfunctions  of 

the  Laplace-Beltrami  operator  on  r,  but  practically  this  sys¬ 
tem  is  difficult  to  construct  explicitly.  Let  us  prove  that  for  a 
starlike  domain  D  the  system  [  Yt  (£ )  | ,  where  £  =  [0,  <*> )  is  a 
point  on  a  unit  sphereS 2  and  J'  are  the  normalized  spherical 
harmonics,  forms  a  basis  in  each  of  H^.  A  domain  is  called 
starlike  if  there  exists  a  point  0  inside  the  domain  such  that 
every  point  of  the  boundary  of  the  domain  can  be  seen  from 
this  point.  This  means  that  the  equation  of  the  boundary  is  of 
the  form  r  =  R  (0,  <f> )  =  R  (£ ),  where  the  origin  is  at  the 
point  0.  It  is  well  known  that 


e-r--L 


rtf  VS' 


4  irr. 


(S’ 


rM) 

2  «  +  r 


J  =  0.1.2,... 


where  r{{.  =  |g  —  g  *|.  The  system  (  Yj],j  =  0,1,2 . forms 

an  orthonormal  basis  of  II  =  L  !(S and  in  any  Hn  (5 J)  the 
scalar  product  in  can  be  defined  as  («,  t/), 

=  [Q  o"  V  Q  0  9»)o.  j(S  l).  and  ( Yn .  Ym  )f 

=  (2n  +  l),(2m  +  lH/v,  Ym)  =  (2n  +  l)’(2m  +  l)^m. 
where  S„m  is  the  Kronecker  delta. 

Lemma  l-.The  system  (  ))  forms  a  Riesz  basis  of 

=  H^F),  provided  that  D  is  starlike,  FC.C~,  and  the  ele¬ 
ments  of  I/q  are  considered  as  functions  of  £eS\ 

Proof.  Consider  the  eigenfunctions  of  the  equation 


L 


(s')ds' 


=  A„  6„ (s),  ru.  =  |s  —  j'|,  seF. 


Jr  4  -ru- 

Since  D  is  starlike  one  can  rewrite  this  equation  as 


(12) 


*jnpjjvr 
4;r|/?  £)-*(£/)! 


(13) 


wherej  =  R  (5  )  is  the  equation  of  the  surface  Fin  thespheri- 
cal  coordinates, £  =  [8,  <j> ),  0, (f )  =  ifi. (R  (f )), ds  =  pj£  )dg, 
Pj[£  )>0,anddg  =  s\n6d9  d(j>.  The  function 
Po\£ )  =  |2?  b  X  R  i  | .  where  X  denotes  the  vector  product  and 
r  =  R(9,t)  is  the  parametric  equation  of  the  surface  F.  The 
system  [  <Py  J  of  the  eigenfunctions  of  the  operator  Q  defined 
in  (13)  forms  an  orthogonal  basis  of  the  weighted  space 
L  J(S J,  /)0(/ )).  Since(*  *)0  <  c,  <p0{f )  <c2  the  normalized  sys¬ 
tem  [  <Pj  pf  1/2 )  forms  an  orthononnal  basis  of  L  2(5 J). 
Therefore  this  system  is  an  image  of  the  system  {  Yt } , 
j  =  0,1,2,...,  under  a  unitary'  transformation  of  L  J(5J):C/J^ 
=  pf  xn0j  or  <Pj  =  p'0nUYj.  The  operator plnU  is  a  bijec- 
tion  of  L 5 (5 2)  onto  itself.  Let  us  introduce  the  operator  J&t 
—  'fijJ  =  0,1,2,...,  where  arc  the  normalized  eigenfunc¬ 
tions  of  Eq.  (12).  The  operator/ defined  on  the  basis  elements 
is  isometric  and  can  be  extended  to  the  isometric  bijection 
J±  J(S5H-£  7[n  =  H.  Therefore  /y  =  Jp'0nUYt,j 
=  0,1,2,..., and  Lemma  1  is  proved  for//  =  For^^Othe 
system  Y{  forms  an  orthogonal  basis  of (5 ?)  and  the  space 
H%  =  H,{F)  is  metrically  equivalent  to  //f  (S1)  because  F is 
C  “  diffeomorphic  to  S 2.  Thus  the  system  |  Yf )  forms  a  basis 
of  //,{/')  for  any  q.  If  instead  of  C“  diffeomorphism  one 
assumes  that  F  is  C '  diffeomorphic  to  S  *,  then  {  Yj  |  forms  a 
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basis  of  H$,  q<,l.  Lemma  1  is  proved. 

(3)  Let  us  consider  another  projection  method  of  solving 
Eq.  (6)  corresponding  to  the  least  squares  method;  namely 


[Ahm  -  v,  A<t>t )  =  0,  1  <j<r, 

(14) 

or 

j,bJmcm=dJ,  1  <j<n, 

M  M  | 

(15) 

where 

bjm={A^.Af>„).  dj  =  [f,  Adj ). 

(16) 

Since  (4^)  is  a  positive  define  matrix  (if  ker/<  =  |0|  which 
we  assume  for  simplicity',  the  system  ( !  5)  is  uniquely  solv- 

able  for  any  n.  This  system  can  be  obtained  from  the  least 
squares  method  as  a  necessary  condition  of  the  minimum  of 
the  functional 

||  Ah,  -  o||3  =  min,  (17) 

or 

|(/+  T]h.  —  tu|J_  ,  =  min,  to  =  /l0~V.  (18) 


Since  /  +  7"  is  a  bijection  of  H_ ,  onto  itself  and  A  0" 1  is  a 
bijection  of//,  onto  //,  _  ,  the  solution  of  (18)  tends  to 
(/+  T)~  'w  as  «— *00  in//_,  ifte//(1,  and  in //,_  ,  ifoe//,.  . 
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Abstract 

Convergence  of  the  T-matrix  scheme  is  proved  under  more  general 
assumptions  than  in  A.G.  Ramm,  J.  Math.  Phys.  23,  1123-1125  (1982), 
and  for  mere  general  boundary  conditions.  Stability  of  the  nume¬ 
rical  scheme  towards  small  perturbations  of  data  and  convergence 
of  the  expansion  coefficients  are  established.  Dependence  of  the 
rate  of  convergence  on  the  choice  of  basis  functions  is  dis¬ 
cussed. 

Dependence  of  the  quality  of  expansions  in  various  spherical 
waves  on  the  shape  of  the  obstacle  is  discussed. 


PACS  03. 80- +r ,  11.20.-e,  02. 30. Jr. 


1 .  Introduction 

1.  Let  V  be  a  bounded  obstacle  with  the  boundary  T.  Consider 


the  following  problem: 

(v^k2)  11=  0  ;  k>0  .  inH, 

(1) 

“lr=  f. 

(2) 

r  ( ~  ika)  0  ;  i — ►  , 

(3) 

where  fi  is  the  exterior  domain,  and  f  is  given.  Later  we  discuss 
other  boundary  conditions  than  (2)/but  the  basic  arguments  and 
conclusions  will  be  similar  to  those  for  problem  (1)— (3) - 

The  corresponding  scattering  problem  is  as  follows:  find  the 
solution  to  Eq.  (1)  satisfying  boundary  condition  (2)  with  f=0 
and  of  the  form  u=uQ+v,  where  v  satisfies  the  radiation  condition 
(3)  and  uQ  is  the  incident  field.  It  is  clear  that  this  problem 
reduces  to  problem  ( 1 ) —  ( 3 )  for  v  with  f=-uQ  on  r.  Therefore,  we 
discuss  in  what  follows  problem  (1)-(3).  There  is  an  extensive 
literature  about  this  problem.  The  existence  and  uniqueness  of 
the  solution  to  this  problem  for  Liapunov  boundaries  are  estab¬ 
lished  long  ago  and  are  available  in  textbodks  now  [l^  .  The  case 
of  nonsmooth  boundaries  was  also  treated  [2] .  Numerical  methods 
for  solving  problem  (1)-{3)  are  known  (finite  differences,  see 
e.g.  [3],  numerical  solution  of  the  boundary  integral  equations 
of  the  second  and  first  kind  [4] ) . 

Our  concern  is  with  the  T-matrix  scheme  [5]  .  This  numerical 
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scheme  was  widely  used  during  the  last  decade  in  the  problems 
of  acoustic,  electromagnetic,  and  elastic  wave  scattering  by 
one  and  many  bodies,  for  scattering  from  periodic  structures 
etc.  [5-1 0]  . 

Nevertheless  the  basic  questions  concerning  convergence  of  the 
scheme,  stability  of  the  numerical  scheme  towards  small  perturba¬ 
tions  of  the  data  remained  open.  In  [ll]  these  questions  were 
discussed  for  the  first  time.  Here  the  results  from  [ll]  are 
strengthened  and  extended. 

2.  Let  us  describe  the  T-matrix  scheme  in  a  general  formula¬ 
tion.  Let  {t^}  be  a  system  of  outgoing  (not  necessarily  spherical) 
waves,  i.e. 

0  \nA,  «> 

rlTF  —  0  ,  r-*  '  (=) 

From  the  Greer: ' s  formula  it  follows  that 


J<VnWdS  ■  ia^dS,  Vn  ,  (6. 

r  r 

where  u  is  the  solution  to  { 1 ) -  ( 3 )  and  N  is  the  exterior  unit  nor¬ 
mal  on  T  (pointing  out,  into  ft) .  Using  boundary  condition  (2) 
one  writes  (6)  as 

J  =  (n  ,  Vn, 

r 

where 

f n  =  If  ciS  y 
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(8) 


9 


The  T-matrix  scheme  consists  of  the  following.  Let  be  a 


linearly  independent  and  complete  system  of  functions  in 


H0=L*(r).  Let 


m 


K'-  t  ci  • 


(10) 


(m) 


where  c^  ,  1<j<m,  are  constant  coefficients,  which  should  be 


defined  from  the  linear  algebraic  system 


m 


c—  r 

pa"jcj  =  f"  •  16ni 

a„j=  i  +„4jds  =  (<$>j,<l'0). 


m 


Oi) 


(12) 


One  obtains  this  system  if  (10)  is  substituted  in  (7)  and  only 
the  first  m  equalities  (7)  are  used. 


3.  Justification  of  the  T-ituitrix  scheme  requires  positive 
answers  to  the  following  questions: 


Q1 .  Is  (11)  solvable  for  sufficiently  large  m? 


Q2. 
Q3 . 
Q4. 


Does  h  ->  h,  Here  h  is  defined  as  in  (9). 

m  H0 


(m) 


Does  •*  c j  ,  m->-«>?  Is  the  convergence  uniform  in  j. 


Does  the  equality  h=  1  c .  $  .  hold,  where  c.  are  defined 

j=1  3  3  3 


in  Q3 ;  and  it  is  assumed  that  the  limits  c^  exist? 


Q5.  How  does  the  rate  of  convergence  depend  on  the  choice  of  the 


systems  {$..},  (^n>? 
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Q 6.  Is  the  numerical  scheme  based  on  the  equation  (11)  stable 
towards  small  perturbations  of  fn  and  the  matrix  a^j? 

Remark  1 .  In  the  literature  [12]  the  following  questions  were 
discussed:  Is  the  set  of  equations  (7)  solvable?  Is  the  solution 
to  (7)  unique?  These  questions  are  easy  to  answer.  The  set  of 
equations  (7)  is  solvable  for  any  system  satisfying  (4)  and 

(5):  take  the  solution  u  to  ( 1 ) — ( 3 )  (which  does  exist)  and  apply 
Green's  formula  to  and  u  to  obtain  (7).  The  solution  to  (7)  is 
unique  iff  the  system  { }•  is  closed  in  Hq=L  (T)  so  that 

=  0  ,  Vn  =>  h-0.  03) 

p 

This  is  equivalent  to  saying  that  any  h£HQ  can  be  approximated 

with  prescribed  accuracy  e  in  the  norm  of  Hn  by  linear  combination 

m(e)  U 

of  the  elements  t  :  1 1  h—  Z  c.  (e)  ||  <e ,  i.e.  the  system  {ij/.  } 

n  j=l  J  3  J 

is  complete.  We  assume  below  that  the  system  is  closed. 

4.  If  one  takes  as  ^  in  (7)  g  (s  ,y )  =g  (s  ,y  ,k)  =  ^^s-y  ~ 

ye-P/  and  does  not  use  equations  (7)  for  n>1,  then  one  gets  the 
integral  equation 


J  cj  (s,ip  his}  dS  -  i  f dS  s  F  (ip  t  tj  eV. 

Actually,  if  (14)  holds  for  yeBcp,  where  B  is  any  ball  lying 
strictly  in  P,  then  (13)  holds  in  V  because  both  sides  in  (14) 
are  solutions  to  the  Helmholtz  equation  in  p  and  therefore^ if 
they  are  identical  in  B  they  are  identical  in  p.  If  one  lets 
y-*'s  ’€  T  one  obtains  from  (14)  the  boundary  integral  equation  of 


the  first  kind 


Ah=.Jq(S,6)h<S)dS 
r J 

where 


=  bis) 


(15) 


b(5)  *  lim  Fl-^  ,  U  £  V  ,  tj  S 


(16) 


If  one  looks  for  a  solution  of  (15)  of  the  form  (10)  and  uses  a 
projection  method  for  finding  cjro^ ,  one  obtains  the  system 


|  A„.  c‘m-1  =  bn  ,  Unsm, 

j=1  J  J 


(17) 


where 


An^  =  (  A  9j  , 

IfjMsJfhdS,  (19) 

and  the  bar  denotes  complex  conjugation.  The  same  questions 

Q1-Q6  can  be  studied  for  system  (17).  In  this  case  nn  plays  the 

role  of  0  ,  but  now  there  is  no  need  to  assume  anything  about  the 
n 

properties  of  nn  in  ft.  In  fact,  ate  defined  only  on  the  surface 
T.  Questions  Ql ,  Q2 ,  Q5 ,  Q6  were  answered  in  [ll]  for  the  system  (17) 
under  the  assumption  that  nn=4>n  and  the  system  forms  a  basis 

of  H  ,.  The  spaces  H  =W?  (D,-«<q<®  are  defined  as  the;spaces  of 
functions  with  q  square  integrable  derivatives  for  q>0  integer,  and 
as  dual  spaces  (spaces  with  negative  norm)  for  q<0.  For  arbitrary 
q<0  they  can  be  defined  as  interpolating  spaces,  or  directly  [13]. 

In  the  present  paper  we  note  that  the  result  and  arguments  in 


[11]  are  valid  under  the  weaker  assumption  that  { <J> ^ }  is  a 
complete  system  of  linearly  independent  functions  (not  necessarily 
a  basis,  see  Sec.  2.3  below). 

From  the  integral  equation  (14)  one  can  go  back  to  the  sys¬ 
tem  (7)  by  assuming  that 


^(Sjlp  S  T.  ;  ^  (20) 

J 

substituting  (20)  into  (14)  and  equating  coefficients  in  front  of 
$ j .  From  this  point  of  view  the  integral  equation  (14)  is  equi¬ 
valent  to  the  system  (7) -(8).  In  the  literature  expansion  (20) 
is  used  with  ip.  being  the  outgoing  spherical  waves  and  being 
the  regular  (i.e.  finite  at  the  origin)  solutions  to  Helmholtz* 
equation.  Matrix  (18)  will  be  identical  to  matrix  (12)  if  ipn 
in  (12)  are  chosen  so  that  Ann=^n*  This  corresponds  to  a  specific 

choice  of  the  outgoing  waves,  since  for  any  linearly  independent 

2  "  - 

system  of  functions  nn  in  L  (T) ,  the  system  ^n=Ann  will  be  a  sys¬ 
tem  of  outgoing  waves  whose  boundary  values  on  T  form  a  linearly 

2  2 
independent  system  in  L  (T)  ,  provided  that  the  operator  A:  L  (T)-*- 

2  2 
-*L  ( r )  has  no  zeros  (i.e.  An=0=»n=0)  .  This  will  be  the  case  iff  k 

is  not  an  eigenvalue  of  the  interior  Dirichlet  Laplacian  in  V. 

In  [ll^j  it  was  noted  that  in  the  case  when  k2  is  an  eigen¬ 
value  of  the  Dirichlet  Laplacian  in  V,  one  can  use,  instead  of 
g(x,y,k),the  Green’s  ‘function  g£(x,y,k)  and  in  this  case  the  cor¬ 
responding  operator  A  will  have  no  zeros.  The  function  g£  is  the 

2  2 

Green's  function  of  the  Dirichlet  operator ,V  +k  in  the  exterior 

2 

of  a  small  ball  B0  situated  in  V,  where  B  is  so  chosen  that  k 
is  not  an  eigenvalue  of  the  problem 
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(v^+k*)  u=  0  in  P\&£  , 

£ 


where  3B  is  the  boundary  of  B  ,  and  tKB  is  the  complement  in 

C  EE 

V  to  B£. 

The  above  argument  shows  that  the  analysis  in  [ll]  is  applic¬ 
able  to  equation  (11)  under  some  special  choice  of  in  (12}. 
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2.  Analysis  of  the  T-matrix  scheme 


In  this  section  we  discuss  questions  Q1-Q6  formulated  in 
Sec .  1.3. 

1.  The  system  (11)  is  solvable  for  a  given  m  iff 

det(a  .  )m  .  .fo.  For  the  following  analysis  we  need  some  defini¬ 
ng  n,  j  =  i 

tions  and  results  from  the  theory  of  Hiibert  spaces.  These  defi¬ 
nitions  and  results  are  given  in  Appendix  1 .  In  Appendix  3  some 
results  about  convergence  and  stability  of  projection  methods  are 
given . 

We  are  interested  in  the  properties  of  the  coordinate  systems 
{<}>..}  and  which  imply  positive  answers  to  questions  01#  Q2 , 

Q3,  and  Q6.  Let  us  write  equation  (7)  as  an  operator  equation 

ah=f  t  a:  H0-*  l  }  H0=  L?(P},  <7*) 

where  the  operator  a  is  bounded  and  defined  on  all  of  iff 


£  Uh,<lof  5  c*  1  lhl2ds  VheH0  ,  c^o. 

n*i  r 

If  this  inequality  does  not  hold  for  all.heHg#  but  the  system 
{^n>  is  a  basis  of  Hq ,  then  a  is  densely  defined  (i,e,  its  domain 
D(a)  is  dense  in  KQ) .  Indeed,  in  this  case  the  biorthogonal 
system  {^n>  is  also  a  basis  of  HQ  [17,  p.  307^,  and  any  linear 


combination  T.  c.  ij».6D(a). 
j  =  1  3  3 

The  operator  a  transforms  a  function  heH^  into  a  sequence 
(h,^n)=/  h  dS ,  1<n<®.  The  range  of  a  is  dense  in  £  provided 


n 


that  for  any  sequence  {d  }€£  the  series  £  d  i|i  (x)  converges  in 

n=1  n  n 

Hq  and  the  system  {ij>n}  is  w-linearly  independent,  i.e. 

CO 

E  d  =0«*d  =0,Vn.  Indeed,  suppose  that  (ah,d) =0 , V  h€Hft ,  where  the 
*  n  n  xi  u 

n— 1  2  oo  _ 

parentheses  denote  the  inner  product  in  l  .  Then  /h  I  d  to  dS=0.  Since 
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oo 

hcHg  can  be  taken  arbitrary  it  follows  that  Z  d^  ^n=0  and 
d  =0#Vn.  The  operator  a-1  will  be  bounded  and  defined  on  all  of 

l2  iff.  || ah j !  >c i  llh!i  /Vh€H0,  c^>0 .  We  use  the  same  notation 

2 

for  the  norms  in  Hq  and  SL  .  This  inequality  can  be  written  as 


EllVMf  >  ct.flhfdS  ,  VK£H„,  c,>o. 

'  r 


n»t 


Therefore  a  and  a  ^  are  both  bounded  iff 


11  h  11 


SUh 

n»t 


iv. 


,  Vhe  H 


0  ) 


ct>  0 


These  inequalities  hold  if f  the  system  forms  a  Riesz 

basis  of  Hq.  Let  us  consider  the  truncated  equations  (7)r  i.e. 

the  system  (11)  as  a  projection  method  for  solving  (7').  Namely, 

let  Q  be  an  orthcsro jection  in  i  defined  by  the  formula 
m 

Qmf =f  ^  =  (f f^  / 0 , 0  ,  and  Pm  be  an  orthoprojection  in  HQ 

on  the  linear  span  of  ($  ^  , .  .  .  <*>m)  •  The  system  (11)  can  be  written 
as 


=  Qmf 


m  i 


m 

s 

j*1 


tm) 


+i 


(IT) 


This  equation  is  of  the  type,  studied  in  Appendix  3  (.see  A3. 2)  . 
Conditions  (A3. 4)  and  (A3. 5)  are  necessary  and  sufficient  for 
equation  (11')  (i.e.  (11))  to  be  uniquely  solvable  for  all  suffi¬ 

ciently  large  n>m0  and  for  the  convergence 


.  llhm  -  h||  0  ,  m  -**>. 


;  (2i) 


Conditions  (A3. 4)  can  be  written  in  our  case  as 
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C  >  O 


(22) 


where  c^...cr  are  arbitrary  constants,  anj  =  (<{>  j  »4>n)  •  Condition 
(22)  can  be  written  as 


cj  Csi  3  "J  3  c  cj  cj’  ^ 

where  one  should  sum  over  repeated  indices  and  the  bar  denotes 
complex  conjugation.  This  last  inequality  means  that  the  following 
matrix  inequality  holds 


(aa")B  \  C  ($)m  ,  Vm>m0)  C>0.  (22,) 

Here  (a)  is  the  truncated  matrix:  (a)  =  (a  . )"*  .  <I>=(4j-»«{>-i)  = 

in  iu  n  j  n  /  j  -  i  j 

=  /<>  -<+>-.  dS  is  the  Gram  matrix  for  the  system  If  A  (4>)  (A  ($)  ) 

p  3  3  3 

denotes  the  minimal  (maximal)  eigenvalue  of  a  self-adjoint  matrix 
0>0,  then  (22^)  holds  if  for  example 


inf  X(  (aa'O  >  >.  ■>  o  sup  A  ( 1  ^  A.  <  00 .  (22„) 

m '  1  m  z 

This  condition  is  convenient  from  a  practical  point  of  view. 

The  conclusion  is  as  follows:  if  ( 2 2 2 )  holds,  then  the  projection 
method  (11)  for  solving  equations  (7)  converges,  i.e.  equations 
(11)  are  uniquely  solvable  for  sufficiently  large  m>mQ  and 
1 1  h m—  h  1 1  -»■  0  as  m-»-tt  ,  where  h  is  the  solution  of  (7).  The  second 
equality  ( 2 2 2 )  holds  for  example  if  the  system  { <f> ^ )  forms  a 
Riesz  basis  of  HQ.  If  we  take  41..=^  on  T ,  then  a  =  4>  and  inequality 
(22.)  holds  iff  the  first  inequality  (222)  holds.  Indeed,  if 


2 

a-t>,  then  (22^)  takes  the  form  (a  m>mg ,  a=a*.  This 

inequality  holds  iff  the  spectrum  of  (a)m  is  bounded  away  from 
zero  by  a  positive  constant.  To  see  this,  let  us  use  the  spectral 
theorem  for  the  self-adjoint  operator  a>0: 


(  {a-  csl  £,40  =  J  Ciz-ci)  d(  , 

0 

,  -  2 

where  E  is  tn=  resolution  of  the  identity  for  a,  min  (t  -ct)=-— 

2  2  t>0  4 

t  ~ct>0  if  t>c.  Therefore  the  operator  a  -ca  will  be  nonnegative 

for  some  c>0  iff  a>ct>0  where  a  is  a  positive  constant  and  in  this 

case  we  can  take  a=c  in  the  inequality  ( 2 2 ^ ) _  Note  that  a=X^ 

where  X  is  the  constant  in  (222) .  If  the  system  { 4> >  is  such  that 

inf X { (0)  ) >a>0  (in  particular,  if  it  is  a  Riesz  basis  of  Hn)  then 


m 


the  conclusion 


as  above  (after  formula  {22^) ) - 


Condition  (A3. 5}  means  that  for  m>m0  the  set  of  vectors 
m 

1  <  j  <-m  is  linearly  independent.  That  is 


m  f 

del  Can^)- >ri..  7  0  (  Vm>m0. 


(23) 


This  condition  follows  from  (22)  (see  Remark  1  in  Appendix  3 
and  formula  (222) ) . 

If  (22)  holds  then  system  (11)  is  uniquely  solvable  for 

all  m>m0,  and  the  function  (10),  where  (c .. ^ }  ,  1<j<m,  is  the 

solution  of  (11),  converges  in  HQ  to  the  solution  of  (7).  The  rat 

of  the  convergence  is  given  by  (A3. 8).  This  rate  depends  on  the 

rate  of  convergence  of  (i-P  )h  to  zero,  i.e.  on  the  rate  of 

m 

approximation  of  the  function  h=a  by  the  linear  combinations 
m 

Ec.<j>..  Stability  of  the  solution  towards  small  perturbations  of 
j  =  13  3 

th6  operator  (i.e.  of  the  matrix  anj)  and  the  right  hand  side  f 


(i.e.  the  sequence  {fn})  follow  from  the  result  1)  in  Section  2 
of  Appendix  3.  Indeed,  consider  the  perturbed  system 


=f-4£"  ■  (,1"> 

J 

where  bni  and  e n  are  the  small  perturbations  of  the  operator  a  and 

the  right  hand  side  f,  respectively.  Let  b^  be  sufficiently  small 

in  the  sense  that  the  operator  b  corresponding  to  the  matrix  b  .  is 

nD 

sufficiently  small  in.  the  norm:  ||h]|  <6.  Here  b:  can  be  con¬ 

sidered  as  an  operator  which  is  defined  as  follows.  The  system 
{^n}  is  perturbed:  I»/n=^n+r|n’  Tkis  perturbation  generates  the 

perturbation  b  of  a  by  the  formula  bh=(h,n  ).  The  matrix  . 

n  nj 

is  then  defined  as  ($j,nn)-  If  |jhj|  < <5  and  6<c  where  c  is  the 

constant  in  (22)  (or  A3. 4),  then  according  to  the  result  1)  in 

Section  2  of  Appendix  3  the  perturbed  system  (11”)  will  be  uniquely 

solvable  for  all  sufficiently  large  m  and  the  corresponding 

m  ,  .  _ _ .  , 

h  =  I  c.  will  tend  to  h=a  f=(a+b)  (f+f  ),  where  f 

m  j  =  1  d  i  €  e  12 

11  f  II  <£•  Thus 


=  (a+bi'j  --  i'f  ♦  la+bT  j£  \\  4c'U*S) 

The  constant  c’  can  be  specified: 


c*  l!(a+b)V  {iflUa'iUKa.bfll.  •  (25) 

Here  we  used  the  identity  (a+b)  ^-a-1  =-  (a+b)  ^ba-1  and  the  esti¬ 
mate  H  (a+b) -a“ 1  ||<||  (a+b)-1  ||  *  |J  a'1  ||  '  ||  b||  ,  which  follows 
-from  the  identity. 

The  estimates  (21),  (24),  (25)  and  the  above  arguments  give 
answers  to  Q1 ,  Q2,  Q5,  and  Q6.  We  now  pass  over  to  a  discussion 
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of  questions  Q3 ,  Q4 . 


(31) 


(see  Appendix  1).  Condition  (31)  implies  also  the  positive 
answer  to  Q4.  Indeed,  if  (31)  holds  then  h  can  be  written  as 


(32) 


and  the  coefficients  c^  are  uniquely  determined  by  the  element 
h.  On  the  ether  hand  we  know  that 


l\hm~Vil\  -*■  0  j  m  -*■  60  (33) 

where 


First, 


we  conclude  that 


because 


(34) 


=  lim  t  s  l,rn  Cj  ^  *  Cj.  <35> 

Secondly,  we  see  from  (32)  and  (35)  that  the  answer  to  Q4  is  yes. 

In  the  above  analysis  the  basic  assumptions  were  (22),  (26), 
and  (31)  and  we  explained  which  of  these  imply  positive  answers 
to^ which  of  the  basic  questions  Q1-Q6. 

2.  In  this  section  we  discuss  the  assumptions  (26)  and  (31) 
and  a  particular  case  of  the  matrix  a^^  for  which  the  convergence 
analysis  is  straightforward.  Note  that  (26)  and  (31  )  deal  only  with 
one  of  the  systems.  Assumption  (22)  deals  with  the  "interaction" 
between  the  systems  and  {tyj}* 

Assumption  (26)  holds  if  the  smallest  eigenvalue  of  the  matrix 

4>i  j  =(<|>i ,  4>  j )  ,  1<ifj<m  is  bounded  away  from  0:  Am>A>0  (see  Appendix 

1).  Assumption  (31)  holds  iff  the  matrix  <1^..  defines  a  bounded  and 

2 

boundedly  invertible  operator  on  i  (see  Appendix  1).  Since 
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this  is  a  self ad joint  matrix,  this  will  be  the  case  iff 


A  £  A  <  *° 

/  v  m  —  J  v  > 


x>o 


where  A  (X  )  is  the  maximal  (minimal)  eigenvalue  of  the  matrix 
m  rc 


1<i/j<^.- 


One  can  measure  the  "interaction"  between  {<{>.}  and  by 

3  3 


the  operator  generated  by  the  matrix  a  .-6  .  =q  ..  The  assumptions 

ng  ng  ng 


S'-Nrvj  »  S'  m  n+.]  , 


<  «S 


(37^  ) 


is  a  compact  operator  on  V 


1  ^  *=° 


(372) 


are  sufficient  for  the  unique  solvability  of  (11)  for  all  suf¬ 


ficiently  large  ,  and  for  the  convergence  in  £*:  |[  c  ^  -c  ||  ■*  0  , 


m*-«>,  where  c  ^  ,0,0. . . ),  c=a  ^f,  a  is  the  operator  on 


m 


£2  with  the  matrix  an j •  Indeed,  assumptions  (37^),  (37^)  and  Fred- 


~  — 1  2 
holm's  alternative  imply  that  a  exists,  is  defined  on  all  of  £ 


and  is  bounded.  This  fact  and  the  special  structure  of  a  imply 
the  above  statement  about  convergence.  To  see  this,  let  us  write 


(11)  as  c  +Pinqc  -f  ,  where  is  the  orthoprojection  in  £2 


onto  the  m-dimensional  space  of  vectors  with  components  c^=0  for 


j>m.  Since  q  is  compact  and  P  -»-I  strongly  in  £  one  concludes  that 


||  q-P  q  ||  ->  0 ,  Therefore  ||  (I+P  q)  ^(I+q-)  1  l|-*0,  m-*«>.  This 


-1 


m 


m  • 


(m) 


proves  the  statement  about  convergence  of  c  to  c. 


2  y 

If  the  system  {$•}  forms  a  basis  of  H_=L  (r.)  and  fe£  , then  the 


CO  2 

solution  of  (7)  is  h~  I  c  .<{> .  where  c  is  the  limit  in  £  when  m->«> 

-i  =  i  3  3 
#m  \  J  ' 

of  the  solutions  c'  ‘  to  (11)  and  the  solution  to  the  equation 


ac=f.  In  this  case  (3?^ )=»(372)  due  to  Fredholm's  alternative  and 
the  uniqueness  of  the  solution  to  the  equation  ac=0.  Let  us  show 
that  ac=0=>c=0.  If  the  system  forms  a  basis  of  Hq,  then  equa¬ 

tions  (7)  and  ac=f  are  equivalent.  But  the  homogeneous  equations 
(7)  have  only  the  trivial  solution  if  the  system  {f  }  is  closed 
in  hq=L2  (T)  .  (This  assumption  about  {t{>n}  is  very  natural  and  was 
made  in  the  very  beginning,  see  Remark  1  in  the  Introduction.) 
Therefore,  ac=0^c=0.  As  our  analysis  shows,  the  behaviour  of  the 
smallest  and  largest  eigenvalue  of  the  matrices  «K,$j)  , 
are  of  basic  importance  in  an  analysis  of  convergence  and  stability 
of  the  T-matrix  scheme.  For  the  operators  of  the  form  a=I+q, 
where  q  is  compact,  the  projection  scheme  is  discussed  below  in 
subsection  4.  In  this  case  the  justification  of  the  projection 
scheme  can  be  easily  obtained.  In  [ll]  the  problem  was  reduced 
to  a  projection  scheme  for  the  equation  of  the  above  form 
(I+T) h=f, where  T  was  compact. 

3.  Let  us  discuss  briefly  the  results  in  [l  l]  from  our 
general  point  of  view.  The  matrix  analogous  to  A^j  in  (18)  in  the 
paper  [l  l]  was  of  the  form  Anj=  /  4>  j)  r  where  the  operator  A 
was  defined  in  (15).  It  was  noted  in  [ll]  that  A=AQ(I+T(k))  where 
Aq>0  in  HQ=L  (F)  and  T(k)  is  compact  in  for  any  -»<q<°°.  Further¬ 
more,  A«:  H  ->H  ,  is  a  continuous  linear  bijection  of  H  onto 

0  q  q+1  4 

H  Therefore,  A  .  =  [(  I+T  (k  ))<}>,  where  [u,v]  s  (A  u,v)  =  (u,v)  ,  and 

q+1  n  j  n  j  u  2 

(u,v)  being  the  inner  product'  in  H  .  The  operator  I+T(k)  was 
q  q 

assumed  to  be  invertible.  This  can  be  done  without  loss  of  gene¬ 
rality  (see  [ll]  and  our  argument  in  the  end  of  section  1.4 
above).  If  Hj)  forms  a  complete  set  of  linearly  independent  func¬ 
tions  in  our  general  argument  shows  that  the  system  (17)  is 

uniquely  solvable  for  sufficiently  large  m,  and  the  answers  to  the 
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remaining  questions  Q2-Q6  are  similar  to  the  ones  given  above 
(see  Sec.  2.4  below). 

In  particular,  we  have  stability  as  m-*-°°  of  the  numerical 
scheme  corresponding  to  system  (17),  with  the  matrix  (A$j,$n) 
and  the  operator  A  defined  in  (15),  with  respect  to  small  pertur¬ 
bations  of  the  matrix  and  {bn>.  In  [l  l}  it  was  assumed  that  the 
system  {<?j>  forms  a  basis  of  This  assumption  is  weakened  here: 

o-nly  completeness  in  H  ,  of  the  system  {<{>.}  is  required.  If  a 
•  “  2  J 

2 

linearly  independent  system  {4>j}  is  complete  in  H^=L  (r),it  will 
be  complete  in  H_^.  Indeed,  Hq  is  dense  in  H_^.  Therefore, for 
any  feH_^  one  can  find  f^eH^  such  that  [[  f -f £  ||  _^<  e  -  If  the  sys¬ 
tem  {*.}  is  complete  in  Hn,  then  one  can  approximate  f  in 

T  *1  ’  *  U  «  n  t  U 

J  m(e) 

by  a  linear  combination:  ||  f  -  I  c  .  (s)  <j>^  ||  n<  e-  Since 

C»  j  “  *|  J  *•  ^ 

II  ull  0- II  u!' -  J  one  concludes  that  ||  f~m^  c  j  (e)  <J>  j  ||  _£<  e .  This  means 

that  the  system  is  dense  in  H_^.  (The  same  argument  shows 

that  this  system  will  be  dense  in  H  for  any  q<0.)  This  remark 

Vi 

simplifies  the  argument  in  [ll]  in  the  case  when  we  do  not  require 
that  the  system.  {?.;}  be  a  basis.  In  [ll]  a  basis  in  L  (T)  was 
constructed  from  the  "distorted  spherical  harmonics”  under  the 
additional  assumption  that  T  is  star-shaped  (i.e.  there  exists 
a  point  in  V  from  which  every  point  of  T  can  be  seen) . 

4.  Let  us  outline  a  proof  of  convergence  of  the  projection 

method  for  solving  the  equation  u+Tu=f  in  a  Hilbert  space  H  under 

the  assumptions  that  T  is  compact  and  (I+T)  **  is  bounded.  The 

projection  scheme  is  as  follows:  the  approximate  solution  u  is 

m  .  . 

sought  in  the  form  u  =  I  cJ  where  { <f> .  )  is  a  complete  set 

j_l  J  J  j 

of  linearly  independent  elements  in  H.  Let  denote  the  ortho- 
projection  onto  the  linear  span  of  { -j  ,  •  -  •  <t>n)  •  The  coeffi¬ 
cients  c!m^  are  to  be  found  from  the  equations 
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(u^+Tu^-f , $ j ) =0 ,  These  equations  can  be  written  as  an 

operator  equation  P  u  +P  Tu  -P  f =0.  But  P  u  =u,  and  therefore 
J  mmmram  m  m  m 

(I+PT) u  =P  f ,  or  (I+T-r_)u  =P  f,  where  r  = (I-P  ) T.  Since  the 
rnmm  m  m  m  mm 

system  { <+> ^  >  is  complete,  ||  P^f-f  || -*■  0  as  for  any  fixed  feH. 

Therefore  I-P  -*-0  strongly.  This. and  the  compactness  of  T  imply 
that  ] | rro | *  —  0 .  Therefore  the  operator  I+T-r^  is  boundedly  inver¬ 
tible  for  sufficiently  large  m: 

ll.T-r„r=  1  U*T)[l-  <l*T)"Vj} 

=  U-Vrj"  V  =  E  (VrJJ  V  i\  «Vrm\\  <  1 . 

j.o 

Here  V=(I+T)  ^ .  One  can  estimate  the  rate  of  convergence  of 

u  =(I+P  T)  P  f  to  u= ( I+T)  ^f.  Indeed 
rnmm 

IU„-  nil «  lUi’P^TAp^-jill  u-T'f'jf  l| 

$  llU*PmT)"lll\?4-fH  -v  lltl-VrjV  -  Vllll-fU 

Let  || V||  <a ,  ||rj|  <em*  and  aem<1,  then 

IHI.W’I  *  S  laeja  -  , 

j=o  m 

ll  U-  VrJV  -  V  11  ^  . 

These  estimates  show  that  the  rate  of  convergence  of  um  to  u  is 
determined  by  e  ,  a,  and  ||  P^f— f  ||  .  The  above  argument  is  a  par¬ 
ticular  case  of  a  known  general  theory  [14] . 

5.  So  far  we  have  discussed  the  case  of  the  Dirichlet  boundary 
condition  (2).  If  one  has  the  Neumann  boundary  condition 

SW  °"r,  <2'' 
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then  our  arguments  are  essentially  the  same:  equations  (6)  lead 

to  equations  (7)  with  h=u|r/  f  =/<»  fdS,  and  now- the  role  of  U; 

1  n  j,  n  '•'n 

in  equation  (7)  is  played  by  the  functions  -^r* 

However,  the  integral  equation  corresponding  to  this  case 
will  differ  from  (15).  Indeed,  in  this  case  from  (2*)  and  the 
formula 


one  obtains 


Jh  =i  gU^fdS  s  Flyl  (  yet)  ,  (14-) 

where  h=u| 

Let  y-*-s'€  r  in  the  above  equation.  Then,  using  the  known 
formula  for  the  limit  value  on  T  of  the  potential  of  double 
layer  one  obtains 

h  =  Bh  -  2b  j  (15-) 

where 

~  f  ~  ^5  q  ( S,  S')  .  „  ~  (  ,  r— 

Bhs  2  ill  (si  “FS;  <3S  ,  bu)  =  lim  rcy) ,  d 6  • ) 

r  -  ^ 

If  one  has  the  impedance  boundary  condition 


-  hik. 
2>N 


(2") 


then  again  from  (6),  one  obtains  (7)  with  h=u|p,  f^/fty^dS  and 


now  the  functions  -  -^p  + 


Play  the  role  of  t|>n  in  equation  (7)  . 


Completeness  of  all  these  systems  in  II=L  (T)  for  the  case  when 
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are  outgoing  spherical  waves, or  any  system  for  which  the 
expansion 


holds, is 


,  lx\<  \y\  , 


easily  seen.  For  example/  if 


J  j^^dS  =  0  ,  Vn, 

r 

then  (38)  and  (39)  imply: 


(38) 


(39) 


1  a  U(s)  ^  dS  =  0  (40) 

r J 

2  2 

Therefore  u (x) =/g (x, s) fdS  solves  the  equation  (V  +k  )u=0  in  V  and 

r 

in  ft,  and  u|p=0.  This  implies  that  u=0  in  ft.  If  u=0  in  ft  and  in 
V  then  f=(-~)  ,  -  {—)  =0.  Here  +  (-)  denote  the  limit  values  on 

(JN  +  ■  o  N  “ 

T  from  the  interior  (exterior) . 


Numerical  experiments 


The  purpose  of  the  numerical  computations  presented  here  is 
to  test  whether  some  commonly  used  complete  set  of  functions, 
e.g.,  outo-rir.g  and  regular  spherical  waves,  also  satisfies  the 
assumptions  made  in  the  previous  sections.  Do  they,  for  instance, 
form  a  Riesz  basis?  Are  they  good  for  expansion  of  functions 
on  r? 

1.  Before  answering  these  questions  it  might  be  illustrative 
to  consider  a  simpler  one-dimensional  case,  e.g.. 


Sinmx  J  ,  *€to/tr]. 


(41) 


This  set  of  functions  is  a  perturbation  of  the  orthonormal  basis 


vf 


!-  sin  mx  by  factors  e  qmx.  Here  the  constant  q  can  be  taken  as  a 

7T 


measure  of  the  eccentricity  of  the  object.  The  motivation  for  our 
choice  of  the  model  example  will  become  clear  from  the  considerations 
given  after  formulas  (45) .  The  Gram  matrix  of  this  model  problem 
can  be  calculated  analytically.  We  have 


-  t<fl  rr  -  r;f  *  {m*n?  1  }  . 


(42) 


VJe  def  ine  the  condition  number  for  an  operator  A  as 

llAll.il-''1  II  • 

Notice  that  the  Gram  matrix  is  always  selfadjoint  nonnega¬ 
tive  and  the  finite  Gram  matrix  is  positive  definite,  provided 
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r  i  i  »m  .«•*  >.  . .  »  •  j-  Kj jrtggjSi ks£i i  A , 


the  functions  { 4>  j }  are  linearly  independent.  Fora  positive  self- 

adjoint  operator  the  norms  I|A||  and  || A- 1  ||  can  be  calculated 

by  the  formulas  ||A||  =A,  ||A-1  ||  =  A-1  ,  where  A=min  t,  A=max  t  , 

tea (A)  tea  (A) 

and  a (A)  is  the  spectrum  of  A.  In  this  case 

K  =  A  A  '  (43) 

Numerical  data  seems  to  show  that  even  after  normalization, 

the  perturbed  system  e  qmx  sin  mx,  m=1,2,3,...,  does  not 

form  a  Riesz  basis  of  L  (  [0  ,  ir] )  .  The  non-normalized  system  { ^ }  is 

not  a  Riesz  basis  because  the  necessary  condition  0<c<inf  [( $  |j 

m>1  111 

for  a  system  {*m>  to  form  a  Riesz  basis  is  violated. 

In  Table  1  we  give  the  condition  number  <  defined  in  equation 

(43)  for  both  {*  }  and  the  normal  ized  functions  {*  /  II*  ||  )  Three 

m  1  m  1 

different  tendencies  are  noticed: 

1)  An  increase  in  condition  number  k  as  the  truncation  size 
grows. 

2)  An  increase  in  condition  number  k  as  the  eccentricity  grows 

3)  The  normalized  functions  have  smaller  condition  number  as 
compared  to  non-normalized  ones 

2.  We  now  consider  the  spherical  waves,  i.e., 

=  h£'Ar)  \  .  ^J) 

’  ,  r  =tx\  . 

Re^nW  =  jjlkr)  Yn(wlJ  '  (44) 

m  \ 

Here  h£  (kr)  is  a  spherical  Hankel  function  of  the  first  kind  and 

j^(kr)  is  a  spherical  Bessel  function.  The  spherical  harmonics 

Yn (u>)  ,  where  w=(0,*)  is  a  unit  vector,  and  8,*  are  the  angular 

2  2  2 

spherical  coordinates,  are  normalized  in  L  (S  ) ,  where  S  is  the 
unit  sphere  (n  is  a  multi-index  n=(£,m)). 


.For  large  orders  (£>>kr)  in  ^  and  it  is  known  [20] 

that 

^n(X)  L(2E-0!!(kr^  1  *Yn 

Re<j>nlX)  -  UzEM)!!) 1  tMfYnUo^.  (45) 

For  large  orders  l  the  spherical  waves  are  essentially  a  pertur¬ 
bation  of  the  basis  (to)  by  a  power  of  kr.  This  is  the  promised 
motivation  for  the  choice  of  the  model  example  in  (41). 

In  the  numerical  data  given  below  we  have  also  included  the 
spherical  harmonics  Y  and  the  functions 

*  Urfl‘'  Yn(.u>\  (46) 

which  are  soli  ns  of  the  equation  Au=0  in  n.  The  factor  k  in 

(46)  was  used  in  order  for  the  factor  kr  to  be  dimens ionle'ss . 

The  Gram  matrix  for  these  four  systems  has  matrix  elements, 
which  are  double  integrals  over  the  unit  sphere.  We  assume  that 

the  equation  of  the  surface  V  can  be  written  as  s=p (w) ,  where 

2  —1 

tue S  ,  and  p  is  a  smooth  invertible  function,  so  that  w=p  (s)  . 

In  bur  calculations  the  functions  (44)  and  (46)  were  considered 
on  T,  i.e.  as  function  of  s,  where  r=  |  s|  ,  x=s  on  the  surface  r. 
If  the  bodies  are  axially  symmetric  then  the  matrix  elements  of 
the  Graraians  ran  be  written  as  single  integrals  in  6,  and  for 
simplicity  we  choose  the  surface  of  a  spheroid,  i.e. 


r(B)  =  (sin2, 9 /a1 


2 

+  COS 


(47) 
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V 


where  a  and  b  are  the  semi-axes  of  the  spheroid  (8  is  the  polar 
angle,  so  that  b  is  the  semi-axis  along  the  axis  of  symmetry).  The 
mirror  symmetry  r(0)=r(ir-0)  implies  that  Jt  even  and  l  odd  do 
not  couple.  Thus  the  elements  of  the  Gram  matrix  are  zeroes  if 
£,  +  £,'  is  odd.  In  this  case  we  can  change  the  enumeration  of  the 
columns  and  rows  in  the  matrix  so  that  it  becomes  a  block  dia¬ 
gonal  matrix  with  two  blocks.  The  size  of  the  first  block,  which 

is 


^ax+1 

corresponds  to  the  rows  with  even  numbers,  is  — = -  if  Z 

j  4  max 

odd,  and  — ~^+1  if  Z  is  even.  The  size  of  the  second  block. 

2  max  .  f 

l - +1 


which  corresponds  to  the  rows  with  odd  numbersfis 

£ 

Z  is  odd,  and  -  if  £,  is  even. 

ITlclX  4  Kiel  X 


max 


if 


The  numerical  computations  seem  to  indicate  that  neither  of 

*  2 

the  systems  {#  } ,  {Re^n>ror  (xn)  forms  a  Riesz  basis  of  HQ=L  (T) . 
However,  {Y^ (p  '  (s) ) }  forms  a  Riesz  basis  of  Hq  as  has  been 
proven  earlier  [ll]  .  It  is  seen  from  Table  2  that  the  spherical 
waves  {Reijjn},  and  (xn)  are  not  good  for  expansions  in  the 

sense  that  the  condition  number  of  the  Gram  matrix  grows  as  the 
truncation  size  increases.  Indeed,  in  this  more  realistic  example 
the  tendencies  1)-3)  discussed  above  for  the  model  problem 
seem  to  be  valid,  together  with  the  additional  observation 


4)  The  systems  {ty  }  and  (x_)  have  smaller  condition 
number  than  {Reij;^}. 

The  numerical  data  seems'  to  indicate  that  the  normalized 
functions  should  be  used  for  expansions  of  the  surface  field 
since  the  corresponding  Gram  matrix  has  a  smaller  condition 
number.  Furthermore,  there  is  an  indication  that  for  high  trun¬ 
cations  the  normalized  systems  {tj^}  and  (xn)  are  better  than 

{ReijJ^}  However,  the  difference  is  not  very  considerable.  A 
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large  condition  number  means  that  the  Gram  matrix  is  difficult  to 
invert  numerically.  It  also  means  numerical  instability,  i.e. 
strong  dependence  of  the  numerical  results  on  the  roundoff  errors 
and  errors  in  the  data. 

It  should  be  noted  that  the  Gram  matrix  is  identical  to  a 

r  -i 

Q-matrix  |_6J  in  the  T-matrix  scheme  for  a  special  choice  of  ex¬ 
pansion  functions.  Thus/ some  of  the  properties  mentioned  above 
might  appear  within  this  scheme  for  this  special  choice  of  the 
expansion  functions.  It  should  also  be  noted  that  when  the 
T-matrix  scheme  is  used  to  compute  the  scattered  field,  additional 
operators  enter,  which  tend  to  improve  the  situation.  However,  a 
discussion  of  these  aspects  requires  further  investigation. 


4.  Concluding  remarks 

There  are  several  questions,  raised  by  the  T-matrix  approach 
to  scattering,  which  require  further  study.  We  conclude  by 
giving  a  short  list  of  such  questions. 

expttklx-sQ 
4tt  ix-sl 

r 

of  outgoing  waves  for  calculations.  Here  the  system  (4>n)  is  a 

2 

complete  system  in  HQ=L  (T)  .  If  the  system  {<f>n>  forms  a  Riesz  ba¬ 
sis  of  then  the  system  {vn>  on  T  forms  a  Riesz  basis  of  HQ. 

This  .was  established  in  [ll].  From  the  results  in  [ll]  it  follows 
that  any  solution  of  the  problem  { 1 )  —  ( 3 )  can  be  represented  by 

the  series  u=  Z  c  v  (x)f which  converges  in  ft  up  to  the  boundary. 

n  =  1  n  n,  2 

Indeed,  assume  (without  loss  of  generality)  that  k  is  not  an 

eigenvalue  of  the  Dirichlet  Laplacian  in  P.  Then  the  solution  to 

the  problem  ( 1 )  —  ( 3 )  can  be  written  as  u  (x)  =  /  "X?  s  ^  o  (s  )dS  , 

r  4lT  lx“®  I 

gxo  ( ik  1  s 1  *"S  1 ) 

where  a  is  the  unique  solution  to  the  equation  /  ■ — 4^  |  s ' -s  } — Li~0(s,)dS= 

=f(s).  If  (<j>  )  forms  a  basis  of  L  iT)  ,  then  a=  T  c  4>  .where  the 
n  _  1  n  n» 

2  00  n-1 

series  converges  in  L  (T).  Therefore  u(x)=  E  c  v  (x), where  the  series 

n=1  n  n 

converges  uniformly  in  the  closure  of  ft.  This  was  the  reason  for 
suggesting  the  system  (vn>  instead  of  the  usual  outgoing  waves  (44)^ 
which  do  not  seem  to  form  a  Riesz  basis  on  non-spherical  surfaces. 

The  other  reason  for  choosing  {vn>  was  that  since  the  expansion 
of  the  solutions  to  Helmholtz'  equation  in  the  exterior  domain  in 
the  functions  vn  converges  up  to  the  boundary,  no  difficulties  with 
the  Rayleigh  hypothesis  arise. 

2)  It  was  noted  in  Section  3.2  that  the  Gram  matrices  correspond 
(for  a  particular  choice  of  expansion  systems)  to  a  Q-matrix.  More 


S  ,  X  £ 


1)  Try  the  system  Vnlx) 


.  I 
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extensive  numerical  experiments  concerning  the  condition  number 


for  various  matrices  a^=  )  (i.e.  Q-matrices)  are  called  for. 

This  would  provide  a  better  basis  for  judging  the  performance  of 
specific  choices  <{k  ,  tyj  and  thus  provide  a  more  detailed  answer 
to  Q5 . 

3)  In  the  present  article  we  have  concentrated  on  "the  null 
field  equations”  [5]  {i.e.  (7),  (11)  etc.)  and  the  question  of 
obtaining  a  solution  on  the  surface  T.  In  the  T-matrix  scheme  one 
furthermore  computes  the  (truncated)  transition  matrix  (T)m  * 
the  form 

mm=  (aVUtij'  us) 

1  2 

where  the  m  x  n matrices  (Q’)m  #  (Q")m  are  similar  to  the  Q  ,  Q 
matrices  in  Appendix  2.  The  exact  (infinite)  T-matrix  is 
independent  of  the  expansion  systems  used  on  T.  However,  the 
approximate  truncated  matrix^ computed  according  to  (48)fdoes  con¬ 
tain  such  a  dependence.  It  is  of  interest  to  investigate  the  rate 
of  convergence  of  truncated  forms  like  (48)  to  the  true,  infinite 
transition  matrix  for  the  scatterer  in  question.  It  is  then  of 
interest  to  exploit  general  constraints  on  the  scattering  matrix 
such  as  unitaritv  and  symmetry  (see  e.g.  W  and  the  contribution 
by  P.C.  Waterman  in  [5]  )  . 

4)  Extend  the  discussion  of  [l and  the  present  work  to  the  case 
of  a  penetrable  scatterer.  Of  particular  relevance  here  is  the 
relation  between  the  convergence  rates  in  the  expansion  used  for 
the  surface  fields  and  their  normal  derivatives  (one  aspect  of 
this  relation  is  treated  in  Appendix  2) . 
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5)  Study  the  convergence  questions  for  the  T-matrix  scheme  for 
scattering  from  obstacles  with  noncompact  (infinite)  boundaries. 

In  this  context  the  work  in  [2l]  will  be  useful.  In  [2l] 
the  scattering  problem  was  formulated  and  solved  for  domains  with 
non-compact  boundaries.  For  the  boundaries,  which  are  locally  Lia¬ 
punov  and  such  that  outside  of  a  sphere  of  arbitrarily  large  but 
fixed  radius  the  points  of  the  boundary  can  be  seen  from  a  point 
located  outside  the  domain  V  (in  this  case  infinite)  ,  it  was  proved 
in  [21]  that  the  Schrodinger  operator  with  decreasing  real-valued 
potential  has  no  positive  discrete  spectrum  and  the  radiation  condi¬ 
tion  selects  a  unique  solution  to  the  Dirichlet  boundary  value  prob¬ 
lem.  The  case  of  the  third  boundary  condition  was  also  treated. 
Furthermore,  the  existence,  uniqueness,  and  properties  of  the  re¬ 
solvent  kernel  G  of  the  Schrodinger  operator  were  studied  in  detail. 
In  particular, global  estimates  of  G(x,y,k)  and  its  derivatives, 
uniform  in  a<k<C,  G<a<C<°°,  were  obtained  for  |x-y|-»-0  and  Jx-y  [-*-«>.  It 
was  proved  that  the  limit  G(x,y,k  +  ie)  as  e-*0 ,  s>0,does  exist  and  is 
attained  uniformly  in  a<k<C.  This  was  done  for  the  boundaries  T  for 
which  p(s,rQ)  (1  +  |s|“)  +  0  as  |sJ->-«,  s^r.  Here  a>n,  where  n  is  the 
dimension  of  the  space,  and  VQ  is  the  boundary  of  the  "canonical” 

domain  (the  boundary  of  a  cone  if  n>2,  and  of  a  wedge  in  the  two- 

ikr 

dimensional  case).  It  was  proved  []2l]  that  G(x,y,k)=~^-  u(v,y,k)* 

*  (1+o(1) )  as  |x|-r*“>,  xr  1=v,  where  the  functions  u(v,y,k)  are  the 

solutions  of  the  scattering  problem  in  the  sense  that  they  solve  the 

Schrodinger  equation  and  vanish  on  T.  Furthermore,  it  was  shown  that 

2 

an  arbitrary  function  feL  .(f))  can  be  expanded  in  a  Fourier  integral 
in  functions  u  and  a  Fourier  series  corresponding  to  the  negative 
discrete  spectrum  of  the  Schrddinger  operator.  If  the  potential  is 
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equal  to  zero  then  the  Fourier  series  part  of  the  expansion  is 
absent.  The  wave  operators  were  constructed  in  [2l]J  with  the 
help  of  the  eigenfunction  expansions. 
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£  subspaces  of  a  Hilbert  space  and  a  condition  for 


invert ibilitv  of  the  mixed  Gram  matrices  ((fr^,^). 

Let  H ^  and  be  (closed)  subspaces  of  a  Hilbert  space  H. 

Then  the  gap  of  H1  and  H2  is  defined  as 

iip,-ptn, 

where  and  P2  are  the  orthoprojections  onto  and  H2,  respec¬ 
tively.  Clearly  0<6<1 ,  0 (H1 #H2) =0 (H2,H1 ) .  It  can  be  proved  Qu]  that 

=  max  |  supllU-R^xll  Sup  l|il-Pt)x|ll 

»*«*!  BKB«f  J  * 

X£  H, 

The  following  facts  hold  ( [l  4]  ,  p.  252-260). 

Lemma  1 .  Let  and  Kn  be  n-dimensional  subspaces  of  a  Hilbert 

space  H  and  6  (G  , H  )<1.  Then  G„  and  H  have  orthogonal  bases 

ii  n  n 

U1 , . . .un  and  V1 . . ,vn# respectiveiy ,  and  (u± ,v^ ,  1<i/j<n, 
where 

VlUl  -  *  \  . 

Lenuna  2 .  Let  {<£j}  snd  {^j},  two  sets  of  linearly  inde¬ 

pendent  elements  of  H.  Let  0<Xm  and  0<um  denote  the  smallest 
eigenvalues  of  the  matrices  ((f^,^)  and  (4>i<f^)#  1<i,j<m,  respec¬ 
tively.  Let  and  Hn  be  linear  spans  of  ^  , - ,4>m  and  i|»1  , . . .  ,  i|i  , 

respectively,  and  A  =(♦.,*.),  1<i»j<m.  Let  0 (G  ,H  )<1.  Then  A 

m  J-  j  mm  m 

is  invertible  and 


where 


ctm  =  {l-  G(Gm>VO) 


2.  Minimal  svsteas. 


A  system  of  elements  is  called  minimal  if  none  of  the  ele¬ 
ments  belongs  to  the  closure  of  the  linear  span  of  the  others. 

A  minimal  system  is  called  strongly  minimal  if 

limA  =A  >0/  where  A  is  the  minimal  eigenvalue  of  the  Gram  matrix 
mm 

m-*00 

Oij^i'V'  1  <i,j<m.  A  system  {^>  is  called  biorthogonal  to 

**• 

the  system  { .  >  if  (<> . ,  <•> . )  =6  .  . .  The  biorthogonal  system  {<{>.} 

3  3  1  xj  3 

is  uniquely  defined  iff  the  system  { ^ }  is  minimal. 

A  system  {$.}  of  linearly  independent  elements  of  a  Hilbert 
space  K  is  called  closed  in  this  space  iff  for  feH  the  condi¬ 
tions  (f  ,<Sj  )  =0,  V  j  ,  imply  that  f  =  0.  .  : 

A  system  { t  ^  of  linearly  independent  elements  of  H  is  called 

complete  in  H  iff  for  any  feH  and  any  given  positive  number  e>0 

m (e)  m (e) 

one  can  find  an  element  Z  c.4>.  such  that  ||  f-  Z  c.$.  ||<e.  Here 

j  =  1  33  j=1  3  3 

the  constants  and  the  number  m (e)  depend  on  e  and  f. 

A  system  can  be  complete  but  not  forming  a  basis  of  H 

(see  n.3  below  for  definitions  of  bases).  Example:  H=L  ([0,l])f 
0<j<».  This  system  is  complete  in  H,  but  is  not  a  basis 
of  H.  Completeness  follows  from  the  Weierstrass  approximation 
theorem.  The  fact  that  the  system  {x-5}  is  not  a  basis  of  H  can  al¬ 
so  be  easily  explained.  Suppose  that  {x-5}  is  a  basis  of  H.  Then 

eo 

(see  n.3  below)  for  any  fell  the  series  f=  Z  c.xJ  converges  in 
,  j=0  3 

L  i[0#l])-  Therefore  f  is  analytic  in  |x|<1  and  cannot  be  an  ar¬ 
bitrary  element  of  II.  In  fact,  from  Muntz’s  theorem  [16]  it  follows 
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that  the  system  { >  is  not  even  minimal  (every  Schauder's 

basis  is  a  minimal  system).  The  Muntz's  theorem  says  that  the 

system  {xP3),  P0=0,  0<p1<p2< .  is  complete  in  L2  ( [o,l]  )  iff  there 

exists  an  infinite  subsequence  p'.  such  that  E  — r  =<*>„ 

3  j  =  1  Pj 


3.  Bases 


A  system  { g  ^  }  is  called  a  Schauder  basis  of  a  Banach  space  X 

if  any  element  xeX  can  be  uniquely  represented  as  x=  jr  c.a.,  where 

j=1  3  3 

the  series  converges  in  the  norm  of  X.  A  system  {<{>_.}  is  called 
a  Riesz  basis  of  a  Hilbert  space  H  iff  <t>j=Thj»  where  {h..  }  is  an 
orthonormal  basis  of  H  and  T  is  a  linear  bounded  and  boundedly 
invertible  operator  (i.e.  T  1  is  bounded  and  defined  on  all  of 


H) .  A  system  {3^}  biorthogonal  to  a  basis  {4^}  of  H  is  also  a 
basis  of  H. 


A  complete  system  {4^}  in  H  is  a  Riesz  basis  of  H  iff  the 

matrix  9  j )  '  generates  a  bounded  and  boundedly 

2 

invertible  operator  on  s,  . 


A  complete  system  { <j> ^ }  in  H  is  a  Riesz  basis  of  H  iff  there 
exist  positive  constants  a1 ,  a2  such  that 


for  any  m  and  any  constants  c ^ ,  1<j<m.  If  {<£..}  is  a  Riesz  basis 

then  sup  lit-  IS  <  a  _ ,  inf  ||4>.j|  >a.  and  similar  inequalities  hold 
3  1<j<°°  3 

for  ||9.||  .  In  particular  sup  ||  4>  .  ||  <°° .  These  results  can  be  found 
3  1<j<»  3 

e.g.  in  [16]  ,  [  1 7]  ,  [19]  .  In  [15]  the  notion  of  the  Riesz  basis 
with  brackets  is  applied  to  some  nonself ad joint  integral  equations 
arising  in  scattering  theory. 
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4.  Tests  for  boundedness  and  compactness  of  a  linear  operator 

on  £2 . 

Let  (a.  . ),  1<i#j<°°#be  a  matrix  which  is  considered  as  a  linear 

13 

2 

operator  A  on  £  .  When  is  A  bounded  and  when  is  A  compact? 


Here  and  below  we  sum  over  the  repeated  indices. 


x  .  a 


UAvxll*=  ^  |ai,aik\  ,Xl<2  Xj"  ^  SuP  M 


Therefore 


-  « 

I! All  £  sup  (  2  S  \  3';',  3i<l] 

Uk*»®  iM  * 


(At . 1 ) 


CO  00 


The  operator  A  is  compact  if  sup  (  £  £  |a  .-a..  |)3<«>  and 

1<k<~  j=1  i=1  13  1K 


•o  W 


sup  72  2  IsnSt-J 

V=v4  J 


0  .  hi  — >  »  . 


(A1 .2) 


Lemma  (Schur) :  Let  a..=a. 


2  2 
h'.l  ->1  . 


.  .=a..  and  sup  E 1  a  .  .  |  <K  .  Then  ||A||  <M, 
^-3  3  ^  i  n  ^3 


Proof:  It  suffices  to  prove  that  1  (Ax,x)  [<Mj|  x[[  .  One  has 

2  .  1 ..  1  2 


|(Ax/x)l<|aijxjx.j<  1  a.  fla.-Mx.l2  <  M  |[xj[ 


2 


5.  Spaces  with  negative  norms  fl3|  . 

Let  H+  and  K  be  Hilbert  space  H+cH  and  h+  is  dense  in  H. 
Let  ueli+,  feH.  Consider  the  completion  H_  of  H  in  the  norm 


where  |j  .  jj+  denotes  the  norm  in  H+  and  (f,u)  denotes  the  inner 
product  in  H.  The  space  H_  is  a  Hilbert  space,  H+cllcH_,  and 
H  is  dense  in  H- . 
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The  basic  equations  (7)  can  be  considered  as  an  interpolation 

problem  in  the  Hilbert  space  H  (H=H  =L2 (T)  in  our  case),  i.e. 

^  / 

the  problem  (h,i»>n)=f  ,  n=1,2,...  . 

Definition  1  :  A  system  {ii^}  is  called  a  Bessel  system  if 
00  2 

£  (  (h,&  )  |  <°°,  whenever  h£H. 
n=1  n 

Definition  2:  A  system  is  called  a  Riesz-Fischer  system  if 

the  problem 

,  "’U.V'i  (A1.3) 

2 

is  solvable  whenever  {f  }el  . 

n 

The  set  of  sequences  {(h,tj;n)},  n>1  ,  heH  is  called  the  moment  space 

M  of  }.  The  cuestions  of  interest  are: 

n 

1)  When  does  a  sequence  {fn>eM?  That  is,  when  is  the  problem 
(h/^n)=fn/  n> 1 ,  solvable? 

2)  Is  the  solution  unique? 

3)  How  to  construct  the  solution? 

The  answer  to  question  1  is  given  by 

Proposition  1.  In  order  that  (A1.3)  be  solvable  and  ||h||  <C  it 
is  necessary  and  sufficient  that 

I 2  a„f„  i  4  C  I!  I  a„<M 

r= » 

for  any  m  and  any  scalars  a^. 

The  answer  to  question  2  was  already  given:  the  solution  of  the 
problem  (A1.3)  is  unique  iff  the  system  { 4>n }  is  closed  in  H.  The 
following  facts  are  useful  [19], 


Proposition  2-  If  (A1.3)  is  solvable  then  it  has  a  unique  solution 
of  minimal  norm. 

Proposition  3.  The  solution  of 

CA1.4) 

of  minima.',  norm  always  exists,  is  unique,  and  can  be  found  by  the 
formula 

[°  ‘V.  —  4m 

K  =  -UeU<k  6  'f1  del  f1  (A1 .5) 

Mm  1  •  m  \  J  : 

Moreover,  if  (A1.3)  is  solvable  and  h  is  its  unique  solution  of 


minimal  norm  the 


Appendix  2.  Discussion  of  some  expansions  occuring  in  the 


T-matrix  scheme  for  a  penetrable  scatterer 


Using  the  Green's  theorem  one  obtains  the  identity 


*n 

.  r 


(A2.1) 


where  +  denotes  the  limit  value  from  the  interior  (see  (44)  for  the 
definition  of  Re^n)  -  Assuming  that  the  following  series  converge 
in  L2(T) 


^  Cn  , 


(A2.2) 


_  y  \  o 

^  ^  > 


(A2.3) 


one  substitutes  the  series  in  (A2.1)  to  obtain 

0  =  II  i  dnJ  Red>r(  6  S  -  c’-  J  <\>  •  ci  s}  Vn 

n‘  r~\  *  n  1  '  ^  » 


(A2.4) 


Truncating  the  series  and  using  the  first  in  equations  in  (A2.4) 

one  obtains  a  linecir  system,  which  can  be  solved  for  or  dj[^. 

n  n 

In  the  case  when  P  and  Reil>  are  defined  as  in  (44),  it  follows 

Tn  n  • 

from  the  formula 


d  5  -- 


(A2.5) 


that  the  truncated  linear  system  can  be  written  as 


_t«> -r(0  /  AWi  vt\-*W  rv 

Q  d  -  CQ  -  XI)  c  -  0 

where  ^  * 1  ^  =  (d!j .  . .  »d^)  ,  c  ^  =  (c1  . .  .c^)  and 


(A2.6) 


T 


^  ,  n>n'  (A2.7) 

r 

(A  is  a  known  scalar,  depending  on  the  normalization  of  i b  and 

n 

Re^n) .  Term  by  term  differentiation  of  (A2.2)  would  imply 

c  =d  ‘ ,  which  is  not  consistent  with  (A2.6).  That  termwise 

differentiation  can  not  be  used  to  obtain  d  (by  taking  d^=c^* 

n  n 

in  (A2. 2) -  (A2. 3) )  is  most  easily  seen  from  the  following  example 
in  which  T  is  a  sphere  of  radius  a  and  in  which  one  keeps  only 
one  term  in  the  series  (A2.2) ,  so  that  the  questions  about  conver¬ 
gence  of  the  series  are  irrelevant.  Take  u  +=hj|1  *  (ka)  Yfi  (w)  .  Since 
(V2+u2)u  =0  for  i  x !  <a  and  ul  _  =h| 1  *  (ka) Y  (u),one  finds  that 

h  '  [Va  \ 

u .  (x)  =-am"  .  -  ;,^(kr)Y  (w)  (we  assume  that  j  (ka)=0).  Therefore 

T  n  ~n 


r,a=  7n^  ^  ^MYn(^|fsa  =  kl^lka^Yn(^ . 


Indeed 


_  hn;ka'l  =  [jnlkafl  [hnlka,jjjt(ka^- Hnlka]jnlka')]  4  0. 


The  numerator  is  the  Wronskian  of  h^  and  j  and  is  not  zero 

n  n  ’ 

The  general  explanation  is  that  even  if  one  can  continue  the 
series  (A2.2)  analytically  inside  V  in  a  neighbourhood  of  T,  this 
continuation  will  be  different  from  the  solution  u  of  the  problem 

(V+k)u.=  C  in  V  J  alp  “  K+  ,  (A2.8) 

so  that  on .  F ,  where  w  denotes  the  analytic  continuation 

of  the  ser:  (A2.2)  in  P.  Indeed,  u  =w  on  T  and  if  ~  on  r 

+  o  N  3N 

then  usw  in  V  by  the  uniqueness  of  the  solution  to  the  Cauchy 
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problem  for  Helmholtz'  equation.  But  this  leads  to  a  contradiction 
since  w  is  singular  at  the  origin  (OeP) . 

If,  instead  of  (A2.2)  and  (A2.3),  we  consider  the  expansions 


n 


n 


(A2.9) 

(A2.10) 


then  the  truncated  matrix  equation  analogous  to  (A2.6)  reads 

uUtl> 

Q  dn  =  Q  c  }  (A2.li) 


where 


ds. 


""  ~  J  ,vJTn  *Ni 

r 


(A2.12) 


(2) 


Thus,  assuming  that  { , } ,  n,n'<m,is  invertible,  we  obtain 


-*12) 

C  =  a 


(A2.13) 


(this  result  is  used  in  the  T-matrix  apparoach  to  scattering  from 
a  permeable  body  Jfi]  ,  [7]  )  . 


A  relevant 

fact 

Lemma  Let 

(Vt>hu  =  0 

m  V 

in  V 

~  4-1  CJ 

y-i  J 

Assume  that 


{<<Alhu.*0,  |*L*o]  =*•  U=0  iJ, 


(so  that  the  Green's  function  G..  for  the  interior  Neumann  problem 


exists  and  is  unique)  and  that  the  {c^  j<m,  have  been  deter- 

mined  so  that  (where  j|  •  |J  is  the  norm  in  Hq=L  (D) 


u  ML  _  MLmv,  ,  c 
11  bH  ^  Z 


(A2.14) 


Then 


Hu-  Umll  <  C£ 


(A2.15) 


Here  and  below  c  denotes  various  constants  depending  on  r. 
Proof :  Let  u-u  =v„  Then  (V^+k2)v  =  0  in  V  and  ||“|l<e.  We  have 


f  r  .  >  ^g,uU,s)~|  ,  • 

vuf8U^tx,s)^  “  v  JdS 


and  also 


V(s^  *  J  G^jCs.s'i  dS  ,  seT. 
r 

The  following  estimate  is  known  [1] 


Ignu,m\  < 


Is  -  s l  > 


With  h=|—  we  then  have 
oN 


(  ffd  Q  S 

|VCS)I  *  c  )  , 

p  '  • 

I 

which  implies 


l\vl\  «C£ 


since  the  operator  S:  ,  (H^=L  (f)) 


(ShV^H  i 


ViisidS 
Is- si 


is  bounded.  In  factfalso  S:  is  bounded.  Therefore 


where  -^g  is  any  tangential  derivative  of  v.  In  the  above  it  was 
essential  that  the  c  ,  1<j<m,were  determined  so  that  (A2.14) 
was  valid.  If,  instead,  the  cj  ' ,  j<m  are  chosen  so  that 


(A2.16) 


and  if  we  assume  in  this  case  that 


|  u.|p  =  0  ]  u=  0  in  Vt 

(so  that  the  Green's  function  GD  for  the  interior  Dirichlet  prob¬ 
lem  exists  and  is  unique)  we  have 
i  if  bShlXiS)  i,  /r% 

s.N0Vixl  =  -  Sn0  J  “iN  dS  >  *  e  17 , 

where  Nfl  is  a  direction  which  conicides  with  the  normal  to  T  on  T 
However,  when  we  let  x  approach  T,  we  do  not  obtain  a  bounded 
operator  on  KQ  in  the  present  case.  In  fact,  the  estimate  (A2.15) 

does  not  imply 


ms  -£f $fas tua  ,  ito-o , 

r 


(A2.17) 


even  in  the  case  u|j,=f£C  .  Here  u  is  the  solution  to  the  problem 


11=  0  ini?  t  u4p=j, 


(A2.18) 


and  <f>j  solve  the  equation 

<|>j  '  0  ini?. 


(A2.19) 


Proof :  If  we  consider  the  function  f£ei/(r),  ||f£||  <e,  then 

m 

II h  -  Z  c-i  ll  <2e  where  h  =f+f_.  Let  u  denote  the  solution  of 


||h  -  Z  c-$.||  <2e  where  h  =f+f  .  Let  u  denote  the  soli 
e  js1  J  .1  c  e  e 

( 1 )  —  ( 3 )  with  f  substituted  by  h  .  One  can  see  that  II' 


Ec^ 


can  be  as  large  as  one  wishes  if  f  is  chosen  appropriately,  in 
3  u 

fact,  can  be  even  not  defined  on  T.  To  see  this  one  can  take 

k=0  and  V  to  be  a  circle  of  radius  1 .  Then 


•o  .  , 

~  ^  r  e 
n««oo 


,  Wn'sJir  i  r<1 

-TT  * 


-  V  i  *n'"l  *in4 

2>N  "  br  “  c~*  Hn£n  r  e 


n=- 

n*o 


CA2.20) 


00 

If  l  ]nhen|=“  the  function  (A2.20)  has  no  limit  in  L2(T)  as 
«=-“  “22 

r-*1-0.  If  E  Inh^^j  =c  then  the  limit  does  exist  and 


H’Tn’II  =2itc  can  fce  as  large  as  one  wants  if  f  is  chosen  approp- 


h=-<» 


riately . 


Appendix  3.  About  projection  methods. 

1.  Convergence  of  projection  methods. 

Let  A  be  a  linear  bounded  and  boundedly  invertible  operator 

from  a  Hilbert  space  H  onto  a  Hilbert  space  G.  Let  P  be  the 

m 

orthoprojection  onto  L^,  where  L^  is  an  m-dimensional  subspace 

of  H'  Lm+  -p '  and  the  sequence  of  the  subspaces  Lm  is  limit 

dense  in  H,  i.e.  for  any  h£H  the  distance  from  h  to  L  qoes  to 

m  3 

zero  as  m+«>.  Let  Q  be  the  orthoprojection  onto  M  ,  where  M 

m  mm 

is  an  m-dimensional  supspace  of  G,  and  the  sequence 

It  is  limit  dense  in  G. 
m 

Consider  the  equation 


Ah=f  (A3 . 1 ) 

and  the  projection  method  of  its  approximate  solution 


QmAPmhm=  Qmf.  (A3. 2) 

The ‘ question  of  when  the  following  statement  is  true  is  then  of 
interest: 

Equation  (A3. 2)  is  uniquely  solvable  for  all  sufficient¬ 
ly  large  m  and  II h  -h||-*0,  m+«. 

m  (A3. 3) 


Here  h^  is  the  solution  of  (A3. 2).  In  the  problem  described  in  the 

rn  /  \ 

Introduction,  h^,=  1  c'  '<)>.,  and  V  is  the  orthoprojection  in 
*“  i  j  J 

2  J "  2 
H=L  (T)  onto  the  linear  span  of  G=£  ,  Q  is  the  ortho- 

»  m  m 

.  2 

projection  in  £  onto  the  linear  span  of  the  first  m  coordinate 
2 

vectors  in  £  ,  i.e.  onto  the  subspace  of  the  vectors  whose  compo¬ 
nents  f  vanish  for  n>m.  The  following  theorem,  which  is  a  parti¬ 
cular  case  of  a  more  general  result  from  [18]  answers  the  above 
question.  ”  A8°  - 


£4 


Vf 


Theorem  1.  (7A3.3)  holds  iff 


HQmAPmK|\  >.cl\PmM\  ,  Vm>m0>  VheH.OO,  (A3. 4) 

and 

,Vm>m0.  (A3. 5) 

Remark  1-  I—  ar.d  Qm  are  projections  onto  m-dimensional  spaces, 

where  m=1#2,-...  (this  is  the  case  we  are  interested  in  in  this  pa 

per)  then  (A3. 4)  implies  (A3. 5),  because  the  operator  Q^AP  : 

P  H-K)  G  is  an  injective  mapping  between  two  m-dimensional  spaces 
mm 

and  therefore  this  mapping  is  surjective. 

Proof ;  1)  (A3-3)  =?>  (A3 . 4-5)  .  If  (A3. 3)  holds  then  (A3. 2)  is 

uniquely  solvable  for  m>m^  and  therefore  (A3. 5)  holds.  Further¬ 
more,  (Q^A?^ )  *" 1  G_  r-*Af  =h,VfeG.  Therefore  ||  (QmAPm)  “1Qm!|  <c<».  Here 
and  below  c  cer.otes  various  positive  constants.  Thus  IJP^hJj  = 

=  |!  (Q  APra)’"1QinQr.A-"!:ihi!  <c  ||QmAPmh(|  ,  i.e.  (A3. 4)  holds..  Note  that 

(0  AP  )_1.  0  AP  =1  where  I  is  the  identity  in  PH  (not  in  all  of 

'  m  xa  rr.  m  m 

H) .  2  )  (A3. 4-5) -->(A3. 3) .  From  (A3. 5)  it  follows  that  (A3. 2) 

is  uniquely  solvable  m>m0.  To  show  that  Jjh^-hU  ->0,  m-*«,  consider 
the  equalities 


QmA[PmW  +  (  I  ~  PnO  n  ^  ~  ^ m f 

which  imply  that 

QmA.Pm(  s  (A3. 6) 

Since  the  sequence  of  the  subspaces  I.  is  limit  dense  in  II  one 


has  (I-P  )h-»0,Vh£H.  Therefore  (A3. 6)  and  (A3. 4)  imply  that 

II  111!  -  II  Pm(Km-PmM!l  $  ell  QmA(I-Pjh|l  —  0  ,  m-»  (A3. 7) 

Thus 

II  h  -  h„M  ^  B  h-  Pmh  SI  -  a  Pm  h  -  hm  11  -0,  m-*0.  (A3.81 

This  completes  the  proof,  which  is  borrowed  from  [18]  (see  also  [22]) 


2.  Stability  of  the  projection  methods. 

Suppose  that  (A3. 3)  holds  for  the  operator  A  in  (A3.1). 

1)  Will  it  hold  for  A+B  where  ||  B 1 1  <5  and  (5>0  is  sufficiently 
small?  The  answer  is  yes. 

2)  Will  it  hold  for  A+3  where  B  is  compact  and  A+B  is  boundedly 
invertible?  The  answer  is  yes. 

The  proofs  can  be  found  in  [l8].  Since  they  are  simple  we  give 

them  here  for  convenience  of  the  reader.  1)  Let  <$=c-£,  i  1 B  R<  5  r 

where  c  is  the  constant  in  (A3. 4),  0<e<c.  Then  ||q  (A+B)P  h||  > 

m  m  - 

>c  Up  h||  -  <$!i  P  h!j  =e||p  h||  .  For  the  case  when  Qm  and  P  are 

—  ra  m  m  mm 

finite  dimensional  projections  onto  m-dimensional  spaces, Theorem  1 
is  applicable.  (See  Remark  1.)  In  the  general  case  it  is  not  dif¬ 
ficult  to  show  that  (A3. 5)  holds,  i.e.  that  the  operator  Qm(A+B)Pm: 

P  H^Q  G  is  invertible: 

m  m 

Qm  ( A  *  B)  ft.  ’  OmAPm[ln-*  (QmAPj'QmBP„], 

and  ji  <QmAF._ >  “  1 0mBPm (|  <  <1.  Therefore  conditions  (A3. 4-5)  are 

satisfied  by  the  operator  A+B  and  (A3. 3)  holds  for  the  operator 

1  —  1 

A+B.  2)  If  3  is  compact  then  ||  (QmAPm)  QmB  “A  Bll  ?  0,  m+°°. 

because  (Q  A?  )  •+  A~^  strongly.  If  A+B  is  invertible  then  so 

m  m  m 

is  I+A^B,  and  ||  P^h+A^BP^  ||  >c1  ||Pmh||  •  Therefore 
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— ■  — -  ■  ■  -**"••  • -  — -  —  -1-  ••  V-  . 


HQm  ( h*i>)  Pmh  II  =  II  Q„KPm  t  P„h  ♦  (Om  APmi'Qro?,P„hl H 


>•  c||  PmK  ‘  (QraAPj’Q0BPmhll  >,  ell  P„li  *  ABPmhll 
-  c II [  ( QmA p„'l’ Qm B  -  11 

>/  ^  II  Pm'n||  f  V  m>  m„. 

Thus,  condition  (A3. 4)  holds  for  A+B.  To  check  condition  (A3. 5) 

one  notes  that  Q  AP  is  invertible.  Q  BP^,  is  compact  and  0  AP  + 

mm  'mm  c  in  m 

+Q  BP  is  one  to  one  by  virtue  of  (A3. 4).  By  Fredholm's  alterna- 
m  m  J 

tive  one  concludes  that  Q^fA+BJP^  is  invertible  and  (A3. 5)  holds. 
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APPENDIX  VI 


I.  INTRODUCTION 

This  note  is  a  continuation  of  Ref.  1,  where  the  follow¬ 
ing  problem  was  considered: 

(_V2-A-J)u  =  0  in  =0.  (1) 

Here  fl  is  an  exterior  domain,  /'  is  its  closed  smooth  bound¬ 
ary,  and  D  =  R3\/2  is  bounded.  Problem  (1)  has  nontrivial 
solutions  if  and  only  if  (=iff)  k  is  a  complex  pole  At,  of  the 
Green’s  function  G  [x,y,k )  of  the  exterior  Dirichlet  problem. 
In  Ref.  1  a  stationary  variational  principle  for  resonances, 
i.e.,  complex  poles  A?,  was  given 

*2  =  st{<Vu,Vu)/<u,u)J,  (2) 

where  st  is  the  symbol  of  stationary  value, 

{u,v)  —  lim  I  expl  —  er\ n  r)B(x)c(x)  dx, 

«— +oJ 

H;  131 

In  Ref.  1  the  test  functions  for  (2;  were  taken  in  the  form 

u.v  =  r~'txp{ikr)  T  Y.rj .  (4) 

J  -  0 

where  n  =  x\x\~\  YJrt  are  the  sphe-isal  harmonics,  cjm  are 
constants,  A  is  a  parameter,  andgf.-:);;0  .s  a  fixed  smooth 
function  vanishing  on  F  and  equal  to  1  outside  of  some  ball 
containing  D.  It  was  not  proved  in  Ref  I  that  the  numerical 
procedure  suggested  there  converges.  The  question  formu¬ 
lated  in  Ref.  1  concerning  ihe  justification  of  the  numerical 
approach  is  still  open.  The  purpose  of  this  note  is  to  formu¬ 
late  another  variational  prv.ciple  for  calculating  the  com¬ 
plex  poles  k9  and  to  prove  the  converge  nee  of  the  numerical 
procedure.  The  method  in  Re".  1  is  similar  to  Ritz’s  method. 
The  method  suggested  in  tb:s  note  is  similar  to  Trefftz’s 
method.  The  advantage  of  this  method  is  that  one  deals  with 
the  compact  operators,  while  in  Ref.  1  the  operator  was  not 
compact.  Our  construction  is  ratural  in  the  framework  of 
the  singularity  and  eigenmode  expansion  methods.2  The 
convergence  of  the  method  will  be  ptoved.  A  result  which  is 
of  general  interest,  as  i:  se  etn>  to  the  author,  is  a  construction, 
of  a  stationary  variation. r  principle  and  a  proof  of  conver¬ 
gence  for  a  class  of  non-  t.'f-  adjoint  s>  nmetric  operators 
[B*  —  B),  which  occur  fr-.quently  in  the  scattering  theory. 


“Supported  by  AFOSR  $002-  A’.ts  subject  classification  47A10,  78A45, 

81F05, 35705. 


II.  A  VARIATIONAL  PRINCIPLE 

The  starting  point  is  the  following  observation:  k  is  a 
complex  pole  of  G  (xjr.A )  iff  the  equation 


Af == J g(s,t,k  )f(t )  c*  =  0,  ImA  <  0, 


gis.t.k )  =  exp(/A  |s  -  t  |)/(4jr|s  -  1 1),  (5) 

has  a  nontrivial  solution.  This  observation  and  some  conse¬ 
quences  are  discussed  in  Ref.  3.  For  the  convenience  of  the 
reader  let  us  note  that 


Jr  <?A, 


where  A',  is  the  unit  outer  normal  to  I"  at  thepoint  f.  If  k  is  a 
complex  pole  of  G  of  order  r  one  can  multiply  (6)  by  (z  —  k )’ 
and  takez— *Aandx  =  seF.  This  yields  Eq.  (5)  (see  Ref.  3,  pp. 
290-29!)  with /#0. 

Let  us  formulate  the  following  variational  principle 
F[f)=a\Af\\  =  min,  ||/||  =  I,  (7) 

where  \f\p  is  the  norm  in  the  Sobolev  space Hp  — 
ll/ll  —  l/|o-  From  the  above  observation  it  follows  that  (7) 
has  solutions  and  the  min  is  zero  if  k  =  A',,  where  A,  are  the 
poles  of  G(xy,  k  ).  If  A'  /A,  theninf11/(,  ,  \Af\x  >0.  Indeed,  if 
there  exists  a  sequence  ||/„||  =  1,  \Af„  | *0.  then 
fn— v/li/H  =  l,  Af  —  0,  and  therefore  A  =  k%  (see  Ref.  3,  p. 
291).  The  only  point  which  is  to  be  explained  is  the  conver¬ 
gence  in  H : /„—»/■  In  Ref.  3  it  is  explained  that  ,4  is  a  pseudo- 
differential  operator  of  order  —  I,  that  is, 

(8) 

Here  a„a2>0  are  some  constants,  —  oo  </><  oo  i  flCC”, 
and  the  fact  that  k  /A,  was  used  essentially:  if  k  y^A,  then 
ker/f  s=(  f-.Af  =  0|  =  jOJ  and  A  maps//,  onto Hp  +  , .  If 
\Af„  ||— »0  and  ||/„||  =  1,  then  (8)  withp  —  I  shows  that 
II  /„  H — *-0.  This  contradicts  the  equation  ||/,(|  =  1.  Therefore 

inf  1,4/1, >0  if  A  f=k  .  (9) 

Consider  a  numerical  method  for  solving  problem  (7). 
Let  { /  |  be  a  basis  of  N, 

f=r^tcj/r  '  n°) 

I 

The  necessary  condition  for  F[f)  to  be  minimat  and 
min  /’(/'"’)  =  0,  |j/'"'||  =  1.  yields: 

Kj<n,  (II) 
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(12) 


where 


«>,  =«>J* \  =  Wm*4ft)t  £lc>lI>a 

/-• 

Thus 


deta„(A)=0  1  (13) 

Let  A  J"  denote  the  roots  of  Eq.  (13).  Our  first  result  is 
Theorem  1:  There  exists  lim„_,  k  Jj*1  =  A,,  and  it,  are 
the  poles  of Green's  function  G  {xjrjc ).  Every  pole  A:,  isalimit 
of  a  sequence  k  jj”,  where  k  are  the  roots  of  ( 1 3).  Conver¬ 
gence  is  uniform  in  q  for  any  finite  interval  1  <q<Q. 

Proof-.  We  will  prove  that:  (i)  Eq.  (13)  has  roots  in  the 
circle  \k  —  A,  |  <  e  for  any  fixed  e  >  0  however  small  if 
n  >  «(*)  is  large  enough,  (ii'  If.*;  >  n{e\  and  there  are  no  points 
A,  in  the  circle  \k  —  z\<e  then  Eq.  (13;  has  no  roots  in  the 
circle  |A  —  z|  <e.  An  important  par;  of  the  proof  is  the  re¬ 
duction  of  the  problem  to  the  problem  with  the  operator 
1  +  T(k),  where  T(k)  is  compact. 

Let  us  fix  e>0  such  that  in  the  circle  [A:(A  —  A,|  <*| 
=  Kt  there  are  no  other  poles.  The  operator /4  =  A  (k  Jcanbe 
written  as 

A[k)=A0[I+T{k)},  (14) 

where 


A0  =  A[ 0),  AS=A0> 0  in  H=.L\r), 


Aaf— 

f  fdt 

Jr  4  vru  ’ 

05) 

2*(A )  = 

A  o  *[A  (A )  —  A0]. 

(16) 

f'-tcjfj.  (22) 

y- 1 

In  particular,  our  assumption  means  that 
(7W+r*  +  r*r,/+f=(/  +  r*)(/  +  7U, 

(/  +  /». f (A- ))/<"'  =  0=>/<"'=0.  I k-k,\<c,  (23) 

where  Pm  denotes  the  projection  in  H  onto  the  fcear  span  of 
I  /•»—»/. )  •  Equation  (23)  says  that  /  +  P„  T(A)  is  invertible 
in  the  circle  |  A  —  A,  ( <  e.  If  n  is  large  enough  this  implies  that 
I  +  T(k  |  is  invertible  in  the  circle  [k  —  A,  ( <e,  because 
(*))|/  +  f  (A)  —  (/ +  P,  f  (A  )||  -^  as  »-►  « .  Tlbis  a  contra¬ 
diction  since  I  4-  /(A,)  is  not  invertible.  Let  us  explain  (•). 
Weneed  to  show  that  ||(/  —  />J7'||—>0asn-»oo.Since7’(A-)is 
compact  it  can  be  written  as  TN  +  BN,  where  |£*||  <dN, 
ds~* 0  as  N-*  oo,  and  TN  is  a  finite-dimensional  operator.  It 
is  sufficient  to  prove  that  ([(/  —  />„)rv(|— *0as»-»oo.  With¬ 
out  loss  of  generality  one  can  assume  that  T„  isaone-dimen- 
sional  operator,  Tfif=  [f,v)u.  Then 

K7 ~  ~  ll(f - 

<11/11 IMI  as  «—  oo,  (24) 

since/*,— *1  strongly.  Thus  the  statement  (i)  is  proved.  Note 
that  the  orthogonality  of  P„  is  not  used  in  (24),  In  order  to 
prove  (ii)  we  suppose  that  for  any  >  0,  e,— *0.Eq.  (1 3)  has  a 

root  k in  the  circle  |A  —  z|  <em  and  show  that  under  this 
assumption  z  has  to  be  a  pole  of  the  Green's  function.  The 
assumption  means  that 

[I  +  7’(A("i)]/<"1  =  0,  l/Mg-1.  k'*-+z.  (25) 


The  operator  A0  is  a  bijection  of  H„  or. to  HP  ,  , ; 

*il/lo<K/l.<*i!  f'K.,bt,b.  —  ct-r.Sl  >  0,  (17) 

while  T[k )  is  compact  as  a  map  (see  Ref.  3  for  de¬ 

tails)  because  A  (A )  —  A0  is  an  operator  with  a  nonsingular 
kernel.  Let  us  rewrite  functional  (7)  as 

/■(/)  -  Mo (/+  T)/\\  --=  tnin,  !(/![  =  1.  (18) 

From  (18>  and  (17)  it  follows  that  the  problem  (7)  is 
equivalent  to 

11/'!  ~  1*  (»*) 

The  matrix  of  the  system  (11)  can  be  written  as 
«*=((/+  T)fmM+T)fJ),  (20) 

where  (...)  denotes  the  scalar  product  which  is  metrically 
equivalent  to  the  scalar  product  in  H.  This  means  that 
<M/>/)o<(/*/)<</2(//)o-  where  <7,  >  0and</2  are  constants, 
feH  is  arbitrary.  In  the  sequel  we  will  not  discriminate  be¬ 
tween  (...)  and  (.,.)o-  This  is  possibie  because  ((/  +  T)f, 

'(/  +  T)f]  and  ((/  +  T)f  !/  +  T  )/)„  attain  their  zero  values 
simultaneously.  The  systir  (It)  can  be  considered  as  the 
system  which  corresponds  .o  the  Ritz  method  for  functional 
*(19)  with  the  test  functior.s  If,  j.  This  completes  the  reduc¬ 
tion  of  the  original  problem  to  tne  problem  with  the  operator 
/  +  T(k ),  where  T(k )  is  a  cc  mpact  analytic-in-A  operator 
function  on  H.  To  prove  (i )  let  u»  assume  that  for  a  fixed  e  >  0 
and  At,  and  all  n  there  are  no  >  oo  is  k  of  Eq.  ( 1 3)  in  the  circle 
l*-M  <  e.  The  system  ( 11 )  with  the  matrix  (20)  says  that 

((/  +  r)/-',(/  +  T)f)  =  0  1  </<n,/<V 0,  (21) 


Since  ||/l',|(|  =  1,  one  can  extract  a  weakly  convergent  in  H 
subsequence  which  is  denoted  again/1"1/1'’1— ■/(— -  means 
weak  convergence).  Since  7"  (A )  is  compact  the  sequence 
T[z)f>  converges  strongly  in  H: 


T(z)f"]—*T(k  )f 

(26) 

On  the  other  hand. 

||r(Aj-r(z)||-*o. 

(27) 

From  (25H27)  it  follows  that 

ll/ll  - 1. 

(28) 

and 

[/+7W=0,  ||/||  =  1. 

(29) 

The  proof  is  complete. 


111.  DISCUSSION 

The  variational  principles  ( 1 9)  and  ( 1 8)  can  be  viewed  as 
the  least  square  method.  Let  us  consider  instead  of  (13)  and 
(20)  the  following  equation: 

det  bjjk )  =  0,1  bJnt^[I  +  T(A  ))/„,/,).  (30) 

Arguments  similar  to  the  ones  given  in  Ref.  3,pp.  192-193 
show  that:  (i)  For  any  e  >  0  and  A,  there  exists  a  root  A  Jf4  of 
Eq.  (30).  such  that  |Af  -  Ai"'|  <e  if  »_>n(e).  (ii)IfA:,;>  is  a 
sequence  of  the  roots  of  Eq.  (30)  and  A  Ir'1-*Af  as«-+oo,  then 
A,  is  a  pole  of  the  Green’s  function.  Equation  (30)  can  be 
viewed  as  a  necessary  condition  for  the  linear  system  of  the 
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Galerkin  method  for  the  equation  (/  4-  T(k  ))/=  0  to  have  a 
nontrivial  solution.  The  Galerkin  equation  is  of  the  form 

(/*"'  4-  T(k  —  0  l<y<«,  (31) 

where  f1’*  is  defined  in  (22).  The  basic  idea  is  that  the  poles  k9 
are  the  points  at  which  the  operator  /  +  T(k )  is  not  inverti¬ 
ble.  These  points  can  be  found  by  the  Galerkin  method,  by 
minimizing  functional  (19)  or  by  some  other  method.  It  is 
interesting  to  note  that  the  Ga’erkin  equation  (31)  can  be 
obtained  also  as  a  necessary  conci  tion  for  the  stationary  vari¬ 
ational  principle 

((/+r(*  ))/,/)  =  st,  ji  /(•  >  0,  (32) 

where  st  means  stationary  •-alae.  This  is  not  true  for  an  arbi¬ 
trary  operator,  but  the  operator  2?  ~i  4-  T{k )  is  a  symmetric 
non-self-adjoint  operator  on  11  =  L  7{T ).  that  is 

~B»=B  or  £|w)  =  «m[?S(u)].  (33) 

Therefore  the  necessary  condition  for  (32),  which  can  be 
written  as 

[Bf,h )  -f  (/?«,/)  =  0  for  all  k=H,  (34) 

yields 

0  -  (Bfh )  +  (  )  =  (Bfh  \  4-  [Bf.k ).  (35) 

Let  h  =  v,  where  veH  is  an  arbitrary  real-valued  function. 
Then  (35)  says  that 

0  =  /?(/+/).  (36) 

Let  /;  =  it;.  Then  (35)  says  that 

0  =  2?  (/-/).  (37) 

From  (36)  and  (37)  it  follows  that  the  equation 

Bf={I+T{k))f=  0,  ji/Si  >  C  (38) 

is  a  necessary  condition  for  (32). 

Our  aim  is  to  show  that  Eq.  (31)  is  a  necessary  condition 
for  the  problem 

(2?/./)=  st.  || /|!  >0.  (39) 

Let  us  take  f  --=  /,n|  and  rewrite  (39)  as 

£  bjmcK,Cj  =  st,  bJm  ~  [Bf„,fj).  (40) 

/.»!  “  I 

In  general  assumption  (33)  does  not  imply  the  equality 
bj„  —  b'nj.  Therefore  the  following  lemma  is  of  use. 

Lemma  I :  Assume  (33)  and 

/j=7p  j  —  1.2,—.  (41) 

Then 


bjm=bmJ.  (42) 

The  proof  is  immediate. 

Proposition  1:  Assume  (33)  and  (41).  Then  a  necessary 
condition  for  (40)  is  the  system  (31). 

Proof-.  The  operator  B  —  I  +  T(k )  satisfies  (33).  From 
this  and  Lemma  1,  Proposition  1  follows. 

Remark  1:  The  results  of  Sec.  Ill  give  a  convergent  nu¬ 
merical  scheme  for  a  stationary  variational  principle  (32) 
with  a  compact  operator  T  satisfying  condition  (33),  i.e., 
symmetric  non-self-adjointness.  Such  operators  occur  fre¬ 
quently  in  the  scattering  theory.  A  simple  example  is  prob¬ 
lem  (1).  There  are  other  examples  in  Ref.  4. 

Remark  2:  A  numerical  scheme  for  calculating  the  re¬ 
sonances  based  on  theorem  1  is  as  follows:  (I)  Calculate  ma¬ 
trix  ajm  by  formula  ( 1 1 ).  (2)  Find  roots  of  Eq.  ( 1 3).  The  corre¬ 
sponding  solutions  of  (1)  can  also  be  calculated  by  this 
numerical  procedure:  Find/*"1  by  formula  (10)  and 
h*"1  =  AfW  is  the  approximate  solution  of  (1),  which  con¬ 
verges  to  the  exact  solution  of  ( 1 )  as  n~ *  oo .  This  exact  solu¬ 
tion  is  of  the  form  u=A  (*,)/,/—  lim/1"*  as  n— *-oo  and  lim 
here  means  the  limit  in  H  —  L  I(/'). 

Remark  3:  For  numerical  calculations  instead  of  princi¬ 
ple  (32)  one  should  use  the  equivalent  principle 

M  (*)/./)  =  st.  H/ll  >0.  (43) 

The  equivalence  of  (43)  and  (32)  follows  from  the  fact  that  the 
necessary  condition  for  (43)  is  the  equation 

A  {k  )/=  A0{I  +  T(k  ))/=  0,  || /||  >  0,  (44) 

which  i$  equivalent  to  the  necessary  condition  (38)  for  (32) 
because  ker  A  „  —  { 0  J .  If  one  takes /=/M  as  in  (22),  then  the 
analog  of  (30)  is 

det{A  (&  —  0  ,(45) 

and  th-  convergence  of  the  numerical  procedure  follows 
from  the  arguments  given  for  Eq.  (30). 
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